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Break on through to the other side’. 
Jim Morrison 


... all prohibitions are made only to be broken, must be broken’. 
A.S. Byatt. Possession 


Введение 


Эти два доклада? посвящены одной очень старой теме и нескольким 
из многочисленных недавних выводов из нее. 


e В первой части мы покажем теоретическую и практическую эффек- 
тивность применения расходящихся рядов на нескольких исторических 
примерах. 


e Вторая часть будет посвящена строгому обоснованию теории расхо- 


дящихся рядов, которым мы сегодня располагаем благодаря совсем 


недавним достижениям“, и их CBA3AM С теорией асимптотических pa3- 


ложений. 


e Мы завершим изложение — в третьей части — несколькими приложени- 
ями к динамическим алгебраическим (или аналитическим) системам. 


1Прорывайся вперед, к другому берегу. 

?...все запреты, созданные только для того, чтобы их разрушить, надо разрушить. 

3Этот текст взят из двух докладов, представленных автором в журналах X-UPS в 1991 году. 
Первая версия появилась в предварительных публикациях Математического Центра Политех- 
нической Школы. 

4C недавних пор эту теорию можно изложить с помощью средств элементарной и клас- 
сической математики. Доказательства зачастую остаются долгими, техническими и тонкими, 
они требуют некоторых новых методов, которые невозможно изложить здесь. 


1. Кое-что из истории расходящихся рядов 


Речь пойдет не о том, чтобы развить полное историческое исследование 
применения расходящихся рядов в математике (эту работу еще предстоит 
сделать), но о том, чтобы на некоторых примерах показать важность пред- 
мета, констатировать эффективность этого применения, выделяя несколько 
методов и оставляя теоретические объяснения для второй части. С этой це- 
лью я решил описать несколько решающих этапов: Лейбниц, Эйлер, Коши, 
Стокс, Стильтьес, Пуанкаре, Борель и Харди. (К нынешним достижениям 
я вернусь позже.) 


1.1. Суммирование расходящихся рядов: на что можно надеяться? 


В работе над этим разделом я использовал замечательный об- 
зор [Ри]. Для начала я процитирую несколько фраз английского математика 
Дж. Э. Литлвуда, взятых из предисловия к посмертному изданию работ его 
друга: 

Заголовок содержит любопытный отзвук прошлого, в том числе про- 
шлого Харди. В 1828 году Абель писал: «Расходящиеся ряды — это изобре- 
тение дьявола, и постыдно основывать на них какое бы то ни было доказа- 
тельство». В последовавший период критического пересмотра результатов 
от расходящихся рядов просто отказались. Затем настало время, когда 
обнаружилось, что с ними все-таки иногда можно иметь дело. Сейчас 
это само собой разумеется, но в начале века эта тема, никоим образом не 
мистическая и вполне строгая, считалась сенсационной, что же касается 
нынешнего заголовка — как бесцветен в нем аромат парадокса и дерзости. 

Г.Х. Харди, Расходящиеся ряды. (H1] 


Чаще всего к сходящимся рядам и их суммам прибегают с целью до- 
казать числовые или функциональные равенства, такие как 
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казать числовые или функциональные равенства, Такие как 


oe 
log2=1 7+ 3 re 
= ee aes ls ate В 
log 2 = о aa Oe 


Интересно было бы упорядочить многочисленные равенства этого типа: 
второе из них, очевидно, лучше первого для приближенного вычисле- 
ния |052, Tak как быстрее сходится. Отсюда стремление обобщить пре- 
образования на случай расходящихся рядов и их возможные «суммы», что- 
бы увеличить арсенал имеющихся тождеств. Именно в таком духе работа- 
ли математики ХУШ века, в частности, Л. Эйлер. Разумеется, чтобы слу- 
жить основанием для таких вычислений, суммирование расходящихся рядов 
должно подчиняться некоторым правилам: грубо говоря, нужно иметь BO3- 
можность заменить — в вычислениях, содержащих обычные операции, — 
ряд на его сумму без противоречий. 

Пусть D — это С-алгебра числовых рядов с комплексными членами 
(сходящихся или нет): есть операции сложения, умножения на скаляр и про- 
изведение по Коши. Обозначим через С подалгебру абсолютно сходящихся 
рядов. Имеем гомоморфизм С-алгебр 


5: С-С, 


CoO 
который сопоставляет сходящемуся ряду о = »` Gy его сумму 5(0). 
n=0 
Мы хотим определить оператор суммирования S* для «каких-то» pac- 


ходящихся рядов. Вот первые правила для 5*, кажущиеся разумными (мож- 
но привести их и в другом варианте): 
(51) Правило регулярности: если 0 сходится, то 


S(O) S(O) то есть S* продолжает 5. 


(32) Правило инвариантности по сдвигу: 


5” ув —=а0-— 5” ры 
п=0 ‘ak 


(s3) S* должен быть С-линейным. 
(54) S* — гомоморфизм по умножению. 
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На практике будет естественным определить приемы суммирования 571, 
применимые к некоторым совокупностям рядов 21: С C Dy С D. Различа- 
ются случаи: 


e 7)! — векторное подпространство D, и 51 удовлетворяет правилам 1, 
2, 3; 


® D, — подалгебра D, и S| удовлетворяет правилам 1, 2, 3, 4. 


Этот последний случай — очевидно, куда более полезный для генерации 
интересных тождеств — наиболее сложен при теоретическом обосновании 


(правило 4 нелегко проверить). site 
Возникает естественная идея заменять числовые ряды >> а» формаль- 
+ оо n=0 
ными рядами >> апт” и затем попытаться «положить 5 = 1». 
n=0 


Тогда мы получаем оператор суммирования 5, определенный Ha диф- 
ференциальной” С-алгебре сходящихся рядов (C{xr}; х?а/ах) и принима- 
ющий значения на дифференциальной С-алгебре (Оу; х?4/4х) ростков 
функций, голоморфных в начале координат. Задача состоит в том, чтобы 
определить дифференциальные алгебры A; формальных расходящихся ря- 
дов: C{x} с Ay C С[2] и операторы 5*: A; —?, удовлетворяющие «ра- 
зумным» правилам: 

(51) Правило регулярности: 51 продолжает 5. 

(52) S¥ — гомоморфизм дифференциальных алгебр. 

(53) Если J — оператор «асимптотического разложения» (раздел 1.7), 
то JS} — тождество на A}. 

Остается решить, где принимает свои значения отображение 


Sj: Ar -?. 


Это не может быть Oo: условие 3 влечет тогда A; = С{т}. В дей- 
ствительности мы увидим, что требуется «поляризация» за счет выбора 
направления 4, исходящего из начала координат. Тогда получим ? = Ад — 
дифференциальной С-алгебре функций, голоморфных на ростках секторов, 
рассеченных направлением d (размыкание и лучи сколь угодно малы). Мы 
увидим, что для фиксированного направления d есть два «естественных» 
оператора суммирования a WS, (боковые суммирования). В то же время 


>C-anre6pa называется дифференциальной, если она снабжена С-линейным эндоморфиз- 
MOM 0, который является производной: 6(ab) = ad(b) + d(a)b. 
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для рядов с действительными членами хочется принять за ? ростки дей- 
ствительных аналитических функций на ростках из В * — {0} или R™ — {0} 
в начале координат. Эта точка зрения сходна CO взглядами Эйлера. Мы уви- 
дим, что она приводит к определенным затруднениям. 

Вернемся к числовым рядам и рассмотрим два примера: 


о. 
a ee ae eee ee | 


Предположим, что эти ряды принадлежат совокупности рядов D,, снабжен- 
ной оператором суммирования S*, удовлетворяющим условиям ($2) и ($3), 
и обозначим суммы наших двух рядов Sg и ©1 соответственно. Имеем So = 


= 1 — So, откуда So = 4 И 


с. = 6a 
ee eee ee 


откуда 251 =1—-1+1—-1+1-1+...4S; = 1. Итак, 


eee hae LS hae sa 


L240 = оО a= 


По-видимому, Лейбниц был первым, кто приписал сумме 1—1+1—1+ 


+1-—1+... значение 5. 


Все время рассматривая наши два примера, можно угадать два способа 
суммирования, кажущиеся разумными и позволяющие оправдать только что 
проведенное чисто формальное рассуждение. Эти два способа — дающие 
возможность больших обобщений — суммирование в среднем и суммирова- 
ние Абеля. 


Суммирование в среднем 


Идея этого суммирования — в попытке построить положительную ме- 
со 


py dé полной массы 1 на подмножествах М и определить сумму ряда >> ап 


n=0 
как 
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] Shatin); 


где з„ обозначает частичную сумму 


п, 
Sn = > ар. 
p=0 


Результат не должен зависеть OT первых членов ряда. Следовательно, мера 
должна обращаться в 0 на любом конечном подмножестве №, что приво- 
дит к противоречию, если предполагать меру о-аддитивной! Итак, прихо- 
дится ослабить требования и заменить понятие меры на понятие илотно- 
сти [FGA]. Если 6 — подобная плотность, обозначим ((Е) массу Е под- 
множества N (если она существует). Проще всего построить такую плот- 
ность, используя среднее арифметическое по Чезаро: распределить единич- 
ную массу между n+ 1 точкой 0, ..., n° и перейти к пределу по п. Средние 


30 tsit+...+S8n 


tn = п ) 


соответствующие ряду 1—1+1—1+1—1+..., стремятся к 5 при п — со 


(Sp, принимает значения | и 0 поочередно). 
Если обозначить через Е множество четных чисел, получим 6(Ё) = 


= lim и»(Е) = 5. 
Ряд 1 -2-+3—4+5-—6+... этим способом не суммируется. Ho, 
применив этот способ дважды, получим в качестве суммы 7 


В более общем случае для каждого целого п, можно распределить еди- 
ничную массу по некоторому конечному подмножеству N и предполагать, 
что если Е конечно, To (FE) = 0. Если, кроме того, предположить, что ии 
сосредоточена на отрезке [1, п|, задача сведется к заданию правильной теп- 
лицевой матрицы, то есть бесконечной треугольной матрицы 


000 0 0 
010 All 0 0 


Qn0 Qn1 але a es Qnn 0 0 


Соответствующая мера равна fn = (бо +... +6n)/(n + 1). 


1. КОЕ-ЧТО ИЗ ИСТОРИИ РАСХОДЯЩИХСЯ РЯДОВ 13 


коэффициенты которой — положительные действительные числа, суммы по 
строкам тождественно равны 1, а столбцы представляют собой последова- 
тельности, сходящиеся к нулю. Тогда матрица-столбец, составленная из ty, 
получается умножением A на матрицу-столбец частичных сумм Sp. 


Композиция двух суммирований в среднем, очевидно, соответствует 
произведению матриц. Например, 


Е в 
Ld 
550 

ЕЕ 
ds 3 3 

соответствует методу Чезаро, a 
100 
3 1 

2 4 4 : 

Coa Vil 2B. 2: 
18 18 18 


— его использованной выше итерации. 
Преобразование Эйлера соответствует теплицевой матрице: 


CO м 5 
sch. SCL 0 0 

T=| : - № 3% 
52? On +. aOR 0 


Это преобразование усредняет ряды путем преобразования 


1 n 
by = onl (Chap Cai CG, | 


Оно усредняет формальные ряды, связанные с гомографическим преобра- 
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зованием (считая 1 неподвижной точкой): 
CO CO 
fla) = Sanz”, Gy) = ony", Gly) = f( —— 
n=0 n=0 1 — 94 


Это преобразование было введено Эйлером как ускоряющее сходи- 
мость. Например, для ряда 


11 
2-3 
OHO дает 
1 1 1 
— — + — ; 
2 222 933 


что ускоряет вычисление log 2. 

Вместо того чтобы использовать меры /и„, с параметром n Е N и конеч- 
ным носителем, можно использовать меры и. с параметром t > 0 и произ- 
вольным носителем. Основным примером служит илотность Бореля, свя- 
занная с семейством мер Пуассона: 


Этот метод, возникший благодаря Борелю (1899), сопоставляет ряду 
CO 


›` Gn в качестве суммы предел (если OH существует) 
n=0 


CO 

Sean. «ju _ ть 

lim e'S(t), ще S(t) = > т 
п=0 

В дальнейшем мы вернемся к этому фундаментальному методу. Применяя 

его к «ряду Лейбница» 1 —1+1—1+..., найдем 


lim е 'cht = i 
t>+oo 2 


Методы Абеля 


Описав методы усреднения, переходим к методам Абеля. Самый про- 
CO 
стой из них основан на теореме Абеля: если ряд >) An сходится, TO его 


n=0 
сумма задается равенством 
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fore) 
lim > ар. 
xl 


x<1 п=0 


Отсюда идея — которую можно найти у Эйлера — для некоторых расходя- 
CO 
щихся рядов принять за определение суммы ряда ›»` а» предел (если он 


n=0 


существует) 


CO 
lim > Oa. 


x<1 п=0 


e Применяя этот метод к ряду Лейбница, находим 


n—0 х<1 


как и ожидалось. 
e Для ряда 1 -2-+3—4+5-—б6-+... находим 


| 


1 —2e¢+327+...= —~+— 
(Lae) 


с пределом i 


Этот метод обобщается следующим образом. Полагаем x = е`", тогда 


x” = е ". Пусть {Аи} — данная строго возрастающая последователь- 


ность строго положительных действительных чисел, стремящаяся к +00. 
CO 

Сопоставляем ряду >) а„ функцию 
п=0 


nm 


a= хз ine А 
n=0 


(если она существует для всех достаточно малых t > 0) и — обязательно 
при условии существования — определяем сумму как 


CO 
lim > арке": 
t—O 


t>0 n=0 
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В дальнейшем мы остановимся на следующих примерах: 


0 если п = 0 
ий м если n > 4 ne) 
0 если п = 0, 1, 2 
a | И п) еслип>2 | Гарди AA tees 


1.2. Функциональное уравнение для С-функции. Ряд Эйлера 


Более детальное рассмотрение изложенных ниже вопросов можно най- 
ти в [Ни [Bar]. Мемуар Л. Эйлера [Eul] начинается так: 
Соотношение, которое я собираюсь здесь рассмотреть, касается 
сумм следующих двух обобщенных бесконечных рядов: 
ee St о и 
1 1 1 1 1 1 1 Г 
Tm Om gm ДР бт т 


Задачей Эйлера было изучение равенства 


о о ae о а и (s — 1)!(2° — 1) 1 
И, cos( Serr), 
о eee Oe (28-1 — 1)лз 2 

ть 


Точнее, он собирался «проверить» это равенство для всех целых 3, для 
8 = 5 из = 5. Для этого нужно суммировать расходящиеся ряды 


ym” — —2™ 4+ 3% а" — 6" 4 7" — BN +... 


Эйлер использует метод Абеля. Для т = Ou | он получил два изученных 


выше ряда с суммами 5 И i соответственно. В более общем случае OH 
получает при k > 1 


ЕО ts 0, 


7 7 ыы —122* —1 
124 1 92k Вов Pre 7 By. 


В случае Res > 0 (-функция Римана определена равенством 


C(s) =1°° 427743 "4.... 
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Сопоставим ей функцию 


Имеем 


и знаменитое функциональное равенство Б. Римана для (-функции запи- 
шется в виде 


(SRS = оО" cos( 5 sm )I(s)n(s). 


Мы узнаем тождество, угаданное Эйлером за сто лет до Римана. 

Заметим, что теоретический подсчет, основанный на использовании 
сумм расходящихся рядов, приводит к строгому доказательству формулы 
в случае целых 3. (Эйлер не намеревался приводить такое доказательство.) 

Метод, которым Эйлер обосновывал возможность вычислений с помо- 
щью расходящихся рядов, во многом взят из его исследования ряда Эйлера: 


CO 


So (-1)"nl. 


n=0 


Это обоснование, по большей части «экспериментальное», очень любопыт- 
но. Много деталей, связанных с этим примером, можно найти в [Bar]. В со- 
временных терминах общую идею Эйлера можно изложить следующим 
образом: пусть дан ряд с общим членом а»; следует найти программу (или 
программы), порождающие Ay. Процитируем Эйлера (1749): 


Но я уже отмечал по другому поводу, что слову «сумма» надо придать 
более широкое значение и понимать под ним дробь или другое аналитиче- 
ское выражение, которое, будучи разложено |в ряд] согласно принципам 
анализа, приводит к тому же самому ряду, чью сумму мы ищем. 

В другом месте он пишет: 

..summa cujusque seriei est valor expressionis illius finitae, ex cujus 
evolutione Ша series oritur’. 

7 


...сумма некоторого ряда — это значение конечного выражения, чье разложение приводит 
к данному ряду. 
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В [Eu2] Эйлер описал четыре (новых) способа для суммирования ряда 
CO 


>’ (—1)"n!. Путем прямого подсчета можно проверить, как согласуются 
n=0 
эти разные приближения. Мы не будем описывать эти способы (что очень 


хорошо сделано в [Bar]). В связи с ними нас больше всего интересует тот 
факт, что соответствующий формальный ряд 


CO 


fle) = У (тан 


n=0 


является формальным решением дифференциального уравнения Эйлера 
м 
У aU =o, 


Среди других методов рассмотрим преобразование в непрерывную сходя- 
щуюся дробь 
ae eee ee к ee 
РТ 2 
и итерация преобразования Эйлера следует из суммирования до наимень- 
шего члена, описанного в следующем параграфе. (Этот метод очень хорошо 
объяснен в трактате по анализу, принадлежащем Лакруа.) Метод, основан- 
ный на дифференциальном уравнении Эйлера, связан с суммированием Бо- 
реля-Лапласа и будет играть основополагающую роль в дальнейшем. (Само 
же суммирование Бореля-Лапласа связано с борелевской плотностью, опи- 
санной выше [Н]].) 
Позднее Харди [Н2]| суммировал ряд Эйлера, используя описанный 
выше метод Абеля: 


оо 
lim (2 О 5 Caylee log п log(log a); 
—0 

t>0 n=s3 


1.3. Эйлер, Коши, Пуанкаре и суммирование до наименьшего члена 


В начале главы УШ второго тома «Новых методов небесной механи- 
ки» [Р] А. Пуанкаре пишет: 


Между геометрами и астрономами есть определенного рода недора- 
зумения в связи со значением слова сходимость. Геометры, озабоченные 
совершенной строгостью и зачастую слишком безразличные к громозд- 
кости запутанных вычислений, с помощью которых они обосновывают 
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возможность, не задумываясь о ее эффективном использовании, говорят, 
что ряд сходится, когда сумма членов стремится к определенному пре- 
делу, даже если первые члены убывают слишком медленно. Астрономы, 
напротив, имеют привычку говорить, что ряд сходится, когда, к примеру, 
двадцать первых членов убывают очень быстро, даже если следующие чле- 
ны должны неопределенно возрастать. Так, в качестве простого примера 
рассмотрим два ряда с общими членами 1000" /п! и n!/1000”. 

Геометры скажут, что первый ряд сходится, причем быстро,...; но 
второй ряд они сочтут расходящимся. . . 

Астрономы, напротив, сочтут первый ряд расходящимся,... а вто- 
рой — сходящимся. 

Оба эти правила корректны: первое — для теоретических исследова- 
ний, второе — для численных приложений. 


Итак, можно говорить о сходящихся рядах «в смысле геометров» или 
«в смысле астрономов». Заметим, что на практике, в приложениях, мож- 
но констатировать, что прежде всего сходящиеся «в смысле астрономов» 
ряды имеют очень быстро возрастающий общий член, хотя изначально он 
убывал. Итак, то, что Пуанкаре рассматривал как возможность, на самом 
деле — правило. Он не отдавал себе отчет в этом явлении — хотя был очень 
хорошо с ним знаком, — приводя определение асимптотических разложе- 
ний. Наибольшие недостатки теории Пуанкаре связаны с этим замечанием. 
Чтобы преодолеть их, надо обратиться к работам других математиков. Уже 

CO 

упомянутый нами ряд Эйлера S~* (—1)"n!2"*! (смежный с рядом, на KOTO- 


п=0 
рый ссылается Пуанкаре) позволит нам начать понимать, что надо делать. 


В конце концов мы получим теорию, точную асимптотическую теорию, 
в которой противоречие между двумя отмеченными Пуанкаре понятиями 
сходимости исчезнет. 

Следующий анализ ряда Эйлера считается классическим (ср., напри- 
мер, Олвера [O]). Пусть 
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Заметим, что 


fi) =a Sle ig ee Sy 


IRn(x)| < ] Пе "тар = тат 1. 
0 


Итак, остаток А„(х) имеет TOT же знак, что и первое из опущенных сла- 
гаемых (—1)"п!т”, и его модуль мажорируется модулем этого слагаемого. 
Рассуждая как в случае знакочередующегося ряда при проверке «особого 
критерия», получаем оценки 


№2р(1) < 1(т) < Рра(т) 


для всех целых р > 1. 

Если x > 0 фиксировано, то это не позволяет (как в случае проверки 
особого критерия для ряда, общий член которого стремится к нулю) по- 
лучить сколь угодно точную оценку для f(x) (здесь общий член быстро 
стремится к +00). Ho мы получаем оценку, качество которой зависит OT х. 
Точнее, мы увидим, что она экспоненциально хороша, когда х «мало». 

Ясно, что при фиксированном т > 0 наилучшее приближение f(x) 
последовательностью f(x) получается (этим методом), когда скачок 


12р+1(%) — fap(x)| = (2p) !a7P _ 


наименьший из возможных, TO есть когда общий член ряда Эйлера 
CO 


`` (-1)"nla"*! — наименьший из возможных по абсолютной величине. 
n=0 


В таком случае есть смысл остановить суммирование на соответствую- 
щем индексе: это и есть суммирование до наименьшего члена. 

Здесь отношение двух последовательных членов равно Nx по абсолют- 
ной величине. Если фиксированное т > 0 «мало», а п меняется, то это 
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соотношение меньше 1 при п < x ~'—1u больше 1 при п > x! — 1. Итак, 
общий член убывает до п = N ~ 4 \, а затем неопределенно возрастает. 
При малых т получаем ряд, сходящийся в смысле астрономов. Качество 


аппроксимации равно Nic’ ~ №М!№М-№. По формуле Стирлинга, 
№! > М2 МММе М, 
и мы видим, что порядок аппроксимации равен 


М2 М№е-М ~ \/2r/xe7V/", 


то есть эксионенциален. 
Следуя концепции, выдвинутой Пуанкаре в 1886 году (раздел 1.7), ска- 
жем, что ряд Эйлера 
CO 
fe) = 2-1)" ne! 


n=0 


это асимптотическое разложение функции 


со t/ 
— ВЫ 
Иа) = | та 
0 


в начале координат. 

На практике — в многочисленных случаях, когда асимптотическое раз- 
ложение имеет действительные знакочередующиеся коэффициенты, — знак 
погрешности совпадает со знаком первого из опущенных слагаемых, а ее 
модуль мажорируется модулем этого слагаемого. Сделанный нами анализ 
можно продолжить. Точно так же можно констатировать, что ошибка столь 
же мала, как экспонента некоторого порядка (то есть как er" при под- 
ходящем целом К). Например, как показал Коши [С], это имеет место для 
ряда Стирлинга 

+00 (—1)"-! В, 
>. о 


r=1 
Большее удивление вызывает TOT факт, что суть предыдущего резуль- 
тата, то есть что суммирование до наименьшего члена обеспечивает экс- 
поненциальную точность (возможно, некоторого порядка К), приложима на 
практике к асимптотическим разложениям с комплексными членами в гораз- 
до более общем случае, чем к знакочередующимся разложениям. Этот факт 
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устанавливался экспериментально в течение нескольких веков, он широко 
используется математиками, физиками и астрономами, хотя и не получил 
удовлетворительного объяснения вплоть до недавнего времени, в результа- 
те систематического изучения асимптотических разложений Жевре (раз- 
дел 2.1). 

На практике установлено, что наименьший член часто получается при 
№ ~ а/х (или N ~ а/т^), где а > 0 зависит от задачи. В задачах, тре- 
бующих тонких оценок, по-видимому, сложно предвидеть точный номер 
наименьшего члена (или наименьшего уже неподходящего члена). Тогда 
берут «псевдонаименьший член» при фиксированном т Е С, выбирая в ка- 
честве Л целую часть 2 (где К и а зависят OT структуры изучаемой 
задачи). 


1.4. Croke и каустики. Феномен Стокса® 


Еще в начале XIX века английский физик Стокс хорошо понимал раз- 
личие между рядами, сходящимися в смысле геометров и в смысле астро- 
номов. О первых он часто говорил, что они «сперва расходятся, а затем 
сходятся», а о вторых — что они «сперва сходятся, а затем расходятся». 
Более того, он заметил тот существенный момент, о котором мы говори- 
ли выше (и который, по-видимому, полностью ускользнул от Пуанкаре): 
суммирование до наименьшего члена в общем случае имеет экспоненци- 
альную точность. И поэтому вычисления с помощью расходящихся рядов 
парадоксальным образом оказываются гораздо точнее и быстрее, чем с ис- 
пользованием сходящихся рядов! Стокс дает очень красивую иллюстрацию 
этого принципа, вычисляя (с помощью расходящихся рядов) каустические 
полосы в волновой оптике. 

Каустики — это огибающие световых лучей в геометрической опти- 
ке. В волновой оптике наблюдатель помещается на малом перпендикуляре 
к каустике, и требуется определить, где находятся полосы, то есть где об- 
ращается в нуль функция интенсивности света. 

Теория — созданная английским астрономом Эйри — позволяет пока- 
зать, что в подходящих единицах измерения интенсивность света Ha попе- 


Можно также говорить явление Стокса. — Прим. ред. 
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речнике пропорциональна квадрату интеграла (интеграла Эйри) 


со 


Ai(z) = т ] cos( 5 + zt) dt, 
0 


где 2 — параметр, зависящий от положения поперечника и обращающийся 
в нуль на каустике. 

В связи с этой задачей ставились замечательные физические опыты: 
Миллер измерил положение 25 первых полос с точностью почти до четы- 
рех десятичных знаков. Следовало сопоставить теоретические величины (то 
есть нули функции Az) с экспериментальными. Первый результат принадле- 
жал Эйри: используя формулы квадратур и таблицы логарифмов с десятью 
десятичными знаками, он получил правильное значение (до четырех деся- 
тичных знаков) для положения первой полосы. Теперь следует заметить, 
что функция Эйри — это решение дифференциального уравнения Эйри 


у’ — 29 =0. 


Этот факт позволяет получить разложение А1(2) в сходящийся ряд в начале 
координат (то есть по возрастающим степеням 25): 


со со 
3n 


Ai(z) = 3 2/8 ar =. и 
n=0 9'nlP (в + 2) n=0 от! (п + 1 


y3ntrl 


Функция Ai — целая, так как область сходимости этого ряда беско- 
нечна. Он сходится в смысле геометров! Используя этот ряд, Эйри находит 
положение второй полосы. Вычисления трудоемки: хотя ряд сходится, но 
сперва он расходится. .. До успеха, достигнутого опытным путем, было еще 
далеко, когда появился Стокс. У последнего была во всех отношениях за- 
мечательная идея: найти разложение функции Ai на бесконечности вместо 
нуля (то есть по возрастанию степеней 1/2). Он получил это разложение 
с помощью дифференциального уравнения Эйри. 

Заметим, что это разложение несколько сложнее, чем просто ряд, и еще 
следует добавить, что оно берется не по Z, а по одной из ветвей 2/4. По- 
следний факт поставил перед Стоксом задачу, на которую ему пришлось 
потратить многие годы. (Решение — открытие явления Стокса — было сде- 
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лано в три часа утра мартовской ночью 1857 года? [Sto3].) Функция Ai(z) 
допускает следующее выражение для своего асимптотического разложения 
на бесконечности: 


1 1 _2 +3 Г(п+ 2\P(n+ 4] 3\n 
= = 
Te de oe ee (—z) a 


Ряд 


n=0 


расходится (в смысле геометров!), HO сходится в смысле астрономов, если = 
не слишком мало. Это позволило Стоксу очень легко вычислить положение 
всех полос с четыремя верными знаками, за исключением первой полосы, 
где было только три верных знака и где Стокс проигрывал в сравнении 
с Эйри (Zz слишком мало). Совпадение с опытом было полным, но тео- 
рия — чрезвычайно эффективная количественно — оставалась теоретически 
необоснованной. 


1.5. Суммирование сходящихся рядов вне их области сходимости: 
Борель, Линделеф, Харди 


ре со 
Пусть f(z) = > anx” — сходящийся ряд с радиусом сходимости 


В > 0. Пусть 50 Е С таково, что |%o| > В. Нам нужно найти сумму 
оо 
nr 
расходящегося ряда >> ашто. 
n=0 
Для этого можно попытаться использовать аналитическое продолже- 
ние суммы f, соответствующей ряду f, вне области сходимости. Это даст 


нам метод суммирования, «хороший» в указанном выше смысле, так как 


Лондон, март, 19/57. Кот из дома — мыши в пляс. Ты кошка, а я бедный мышонок. Я делал 
то, что, как я думаю, ты бы не позволила мне сделать во времена семейной жизни, я сидел до 
трех часов утра, мужественно сражаясь с математическими сложностями. Несколько лет назад 
я атаковал интеграл Эйри и после сурового разбирательства свел его к легко вычисляемому 
виду. Но там была сложность, которую я не смог преодолеть, хотя почти заболел от попыток 
сделать это, и наконец мне пришлось забросить ее и признать себя неспособным ее победить. 
Несколько дней назад я вновь принялся за нее, и после двух- или трехдневного сражения, 
причем в последний день я сидел до трех, я наконец победил ее... 
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аналитическое продолжение сохраняет суммы, произведения и производ- 
ные и является инъекцией. 

Пытаясь сделать этот метод более точным, мы немедленно сталкива- 
емся с проблемами существования (возможно, не существует непрерывного 
пути ^/, связывающего начало координат с Xo и такого, чтобы функцию f 
можно было аналитически продолжить вдоль Y; в худшем случае f вообще 
нельзя продолжить вне ее области сходимости) и единственности (аналити- 
ческие продолжения вдоль разных путей могут дать разные результаты). 

В этом разделе мы ограничимся только продолжениями вдоль отрезка 
y = |0, xo], связывающего начало координат с хо. Используя эти продолже- 
ния, продолжим функцию f аналитически на наибольшую (включающую 
область сходимости) открытую звезду (относительно начала координат), 
называемую звездой Миттаг-Леффлера для ряда f. Обозначим эту об- 
ласть Et( fay Обозначим (f, Et( f)) аналитическое продолжение f (опре- 
деленной на области {|| < R}) на Et(f). Сразу проверяем — то есть 
предоставляем читателю возможность проверить, — что оператор «сумми- 
рования» 


fr (1, Et(f)) 


— инъективный гомоморфизм дифференциальных алгебр. 


Все вышесказанное снабжает нас (в некоторых случаях) хорошим те- 
+00 

оретическим методом суммирования ряда >’ аито, но, очевидно, хоте- 
п=0 

лось бы располагать «явным» методом вычисления соответствующей сум- 


мы f(%o9) (то есть одним или несколькими алгоритмами ее вычисления, 
лучше всего — с возможностью программирования на машине и приво- 
дящими к результату, точному в разумных пределах, за разумное время). 
Я изложу некоторые из этих методов, не занимаясь здесь их вычислитель- 
ной эффективностью. 

Сперва заметим, что само аналитическое продолжение (если вспомнить 
его определение) приводит к явному численному алгоритму, который можно 


использовать на ЭВМ [СС]. Алгоритм сводится к конечному числу этапов 
== оо 

следующего типа. Ряд >> а„(т — 11)" сходится с радиусом сходимости 
n=0 

В1 > 0. Возьмем тогда т, где 12 — в. < R,. В окрестности х2 имеем 


Ye oe) ->ь (х — 12)", 
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и остается вычислить OD, как функции от а». Ясно, что это возможно, но 
при этом используются сходящиеся ряды. Так что получается ступенчатый 
процесс. 

Перейдем к другим методам. Кажется разумным использовать абелевы 
методы суммирования, о которых мы говорили выше. Пусть A = {A,,} — 


данная строго возрастающая последовательность строго положительных 
oe) 


действительных чисел, стремящаяся к +00. Ряду »` а» сопоставим функ- 


n=0 
ЦИЮ 


Tg) = хз ane on 
n=0 


(если она существует при достаточно малых t). Теперь остается изначаль- 
но выбрать последовательность {Ли } так, чтобы ряд, определяющий f;(x), 


a 00 
сходился. Зная о сходимости данного ряда f(x) = >) anx”, получаем нера- 
п=0 
венства вида 
lan| < САП 


для подходящих C' > 0 u A > 0. Сразу находим простую последователь- 
ность Л такую, что ряд, определяющий }+(10), сходится для всех t > 0 
и всех 20 Е С: 


ое если п > 0 
Ge 
0 иначе 


| (Линделеф). 


У этого метода есть варианты: 


Oo n 
: AnXo 
lim 


ear 2. Г(1 + 6) (Миттаг-Леффлер) . 


li Ds Ee) age (Леруа). 


с Г + п) 


п=0 


Имеем следующие результаты: 


Теорема (см. [H1], 4.11, р.77-79). Пусть р = ›. Gn — сходящийся 


n=0 
ряд. Пусть то — точка, принадлежащая звезде Миттаг-Леффлера для f. 
Обозначим f (xo) значение аналитического продолжения (обычной) суммы 


ряда } в этой точке. Тогда 
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(i) Если Л. = nlogn при п > 0 и Ло = 0, тогда предел 


CO 
lim > м 
t—0r 

n=0 


существует и равен f (xo). 
(ii) Предел 


существует и равен }(то). 
(iit) Предел 


существует и равен }(то). 


Иначе говоря, методы Линделефа, Миттаг-Леффлера и Леруа позволя- 
ют суммировать сходящийся ряд f в его звезде Миттаг-Леффлера. 
Ясно, что эти методы применимы при как максимум геометрическом 


+00 

возрастании а». Они неприложимы к ряду Эйлера У` (—1)” п!” 1. Пы- 
n=0 

таясь суммировать расходящийся ряд этого типа методами Абеля, нужно 


выбирать последовательность Л, возрастающую «гораздо быстрее». Следу- 
ющая последовательность (подобранная Харди [H2], 1941) хорошо приспо- 
соблена к ряду Эйлера и рядам с аналогичным типом возрастания коэффи- 
циентов !° (ряды Жевре (Gevrey), раздел 2.1): 


— nlogn(log(logn)) если п > 2 
"10 если n = 0, 1, 2. 


Прежде чем проверять эффективность соответствующего метода для 
расходящихся рядов, естественно будет проверить его на сходящихся рядах! 
Вот что позволяет утверждать следующий результат (Харди): 


a со 
Теорема (см. [H2]). Пусть f = У». an — сходящийся ряд. Пусть то — 

п=0 
точка, принадлежащая звезде Миттаг-Леффлера для f. Обозначим f (xo) 


1 Возрастание вида |an| < С(п!)° Ап, где С, A, s > 0. 


28 РАСХОДЯЩИЕСЯ РЯДЫ И АСИМПТОТИЧЕСКИЕ ТЕОРИИ 


значение аналитического продолжения (обычной) суммы ряда } в этой 
точке. Выберем Л„ = nlognilog(logn), если п > 2, и An = 0, если 
п = 0, 1, 2. Тогда предел 


CO 
lim > On Tee tan 
t—O 2 0 


t>0 п=0 


существует и равен }(то). 


В дальнейшем мы увидим, что этот метод суммирования Харди (недав- 
но вновь открытый и обобщенный профессором Жюрка (Jurkat)) очень эф- 
фективен: с его помощью находятся суммы «мультисуммируемых» рядов 
в их «обобщенной звезде Миттаг-Леффлера» [J]. 

Предыдущие две теоремы сперва были доказаны в случае геометриче- 
ской прогрессии f(x) =1l4+2+27+...c¢ суммой (1 — 1)". Затем можно 
легко перейти к общему случаю, используя теорему Коши. 

В дальнейшем мы увидим, что многие решения аналитических функ- 
циональных уравнений в формальных рядах мультисуммируемы (например, 
в случае решений аналитических дифференциальных уравнений). К сожа- 
лению, формальные ряды, являющиеся решениями аналитических 4-диффе- 
ренциальных уравнений (даже линейных), имеют вид не Жевре, но только 
9-Жевре (раздел 3.4). (Имеются оценки вида |а»| < Са" Ап са> Ти 
А > 0.) Итак, эти ряды не мультисуммируемы. Тем более ясно, что к ним 
нельзя применить метод суммирования Харли. Тогда можно попытаться вы- 
брать последовательность A, возрастающую еще быстрее, например, A, = 
= ип? или, в более общем виде, A, = ип? (си > Ou a > 1). Соответству- 
ющие абелевы методы суммирования недавно были изучены А. Фрушаром 
(А. Fruchard) [F2]. 

К сожалению (как предостерегал Харди по случаю абелевых методов 
суммирования, определенных для слишком быстро возрастающих последо- 
вательностей Л), соответствующие методы суммирования уже не позволяют 
находить сумму сходящегося ряда в его звезде Миттаг-Леффлера, но, в об- 
щем случае, только в гораздо меньшей области. Это видно уже на примере 
геометрической прогрессии: получается не вся звезда (здесь: С — |1; +00]), 
но только область, в которую переходит левая часть полосы {0 < Imz < 27} 
без двух прямых Ве""/24 из + Re’™/2% под действием экспоненты. (Можно 
получить и всю звезду Миттаг-Леффлера, устремляя параметр а к 1 при 
фиксированном то, принадлежащем звезде.) 
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Вспомним теперь о методе суммирования по борелевской плотно- 
сти, изложенном в разделе 1.1, в применении к целому сходящемуся ряду 


о +00 
f(x) = >». апт" с радиусом сходимости В > 0. Пусть т — комплексное 


п=0 
число, удовлетворяющее условию |x| < R. Для $ > 0 положим 


n оо 
= = te en 
9» = 2, 202” S(t) = 2. 15» F(t) =e "54. 
p= n= 


Имеем В 
= n+1t” 
Fie (SG) =S@)) Se" У ата ee 5 
n=0 
откуда 


| 
Q 
i=) 
+ 
Oo 
. 
is 
м} 
8 
Q 
3 
+ 
= 
3 
= 
| 


| 
Q 
=) 
+ 
==. 
© 
> 
8 
Q 
= 
+ 
Qu 
In, 


Отсюда легко выводим, что 


+ со м +00 n 
2 ant” = _ jim F(t) =ao + | е—&/= т d& 
n= 0 n= 


Чтобы упростить дальнейшее изложение, предположим, что Ag = 0 
(отсюда легко перейти к общему случаю). Принимая во внимание предыду- 
щие вычисления, приходим к рассмотрению борелевского преобразования 
формального ряда f: B f= $, где 


+00 en-l 
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Ряд ф имеет бесконечный радиус сходимости; поэтому его сумма ф — целая 
функция. Можно проверить, что скорость ее возрастания экспоненциальна, 
не более чем первого порядка. Имеем 


+00 i 
F(t) = Sana” = СФ) = || ети a. 


Итак, можно найти сумму f ряда в его области сходимости как преобра- 
зование Лапласа целой функции 0. 

В более общем случае можно заменить «контур интегрирования» R* = 
= |0, +oo] в интеграле Лапласа на произвольный луч 4, исходящий из 
начала координат. Получаем 


Сейчас мы увидим, что этот формализм часто можно использовать, 
получая также сумму ряда fB некоторой части его звезды Миттаг-Леффле- 
ра, вне области сходимости. 

Для т € С назовем областью Бореля, соответствующей точке т, за- 
мкнутый круг с диаметром [0, x]. Обозначим его О. Получаем следующий 
результат: 


Теорема ([Во1]). //Гсть = сходящийся ряд. Если то — отличная от 
нуля точка звезды Миттаг-Леффлера ряда f такая, что область Боре- 
ля О.о, содержится в звезде ЕК }), и если @ — сумма формального боре- 


левского преобразования ряда f, а а— луч, исходящий из начала координат 
и содержащий то, то интеграл 


La(d)(x) = ] (еде = de 
а 


сходится и равен f (x9). 


К сожалению, процедура суммирования методом Бореля-Лапласа не 
позволяет в общем случае вычислить f (xX) в каждой точке то звезды Мит- 


таг-Леффлера ряда f, но только в подходящей открытой выпуклой области, 
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лежащей в звезде и содержащей область сходимости. В простейших случаях 
получаем выпуклый многоугольник (ограниченный или нет). Например, для 
геометрической прогрессии находим полуплоскость {Rex < 1}. 

На самом деле ряд f легко суммировать во всей его звезде Миттаг-Леф- 
флера, если в методе суммирования Бореля-Лапласа ввести переменный 
параметр и получить так называемый метод суммирования Бореля-Лапла- 
са К-го уровня. Пусть К > 0. Введем оператор разветвления, рассматривая 
функции (определенные на римановой поверхности логарифма) рк(})(5х) = 
= f(a'/*). Положим B. = (pr) Brpe и Ска = (pe) Lax Pk (где направ- 
ление d* соответствует d для разветвления pz). 

Обозначим Г) образ борелевской области D при конформном преоб- 


разовании x => x!/*, Если 0 < k < 5, этот образ берется на римановой 
поверхности логарифма. При k = 2 образ D2 ограничен полулемнискатой 
Бернулли. Мы будем говорить, что D; — это А-область Бореля, и А-область 


Бореля «диаметра» [0, хо] обозначим Ок хо. 


Теорема. /Густь k > 0, a fos сходящийся ряд. Если то — отличная 
от нуля точка звезды Миттаг-Леффлера ряда f такая, что К-область 
Бореля О, содержится в звезде Et(f), и если @ — сумма формального 


борелевского преобразования Ву (1) ряда f, а а — луч, исходящий из начала 
координат и содержащий то, то интеграл 


Ск;а(Ф) (x) 
сходится и равен f (xo) 


При 5 < К угол К-области Бореля в начале координат равен л/К. Чем 


больше К, тем сильнее «разрежены» К-области Бореля. Легко видеть, что, 
так как то — фиксированная точка звезды Et(f), всегда существует действи- 
тельное kg > 0 такое, что для всех k > Ко К-область Бореля содержится 
в звезде. Тогда можно вычислять (то), используя интегральный метод Бо- 
реля-Лапласа произвольного уровня К > ko. 

Если обозначить 5 оператор суммирования сходящегося ряда в его 
области сходимости, то только что описанные нами методы суммирования 
выражаются через операторы суммирования (по направлению (]): 


ба = Са5В (Борель-Лаплас); 
Sed Lege В, (Борель-Лаплас уровня К). 
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1.6. Борель и Стильтьес 


В 1886 году в [Stil] Стильтьес изучает суммирование ряда Эйлера 
при отрицательных значениях переменной. Это куда сложнее, чем в случае 
положительных переменных, так что теперь речь идет о суммировании рас- 
ходящегося ряда, все члены которого положительны. Стильтьес предлагает 
следующий результат: 


1 1 2h: ks В = 
о 


а 
где e* = 1010 (то есть а <> 23, 025851), и показывает, что качество аппрок- 
симации экспоненциально порядка e~%\/27/a, то есть 10710. 

Достоинство метода Стильтьеса в том, что он наделяет ряд из действи- 
тельных членов действительной суммой. Для метода суммирования Бореля 
это не так: для него отрицательное направление действительной оси син- 
гулярно, и там есть два разных оператора суммирования St, и S p-. 
Можно проверить, что сумма Стильтьеса — не что иное, как среднее ариф- 
метическое двух сумм Бореля. Суммы Бореля отличаются от суммы Стиль- 
тьеса на мнимое число порядка 107 105. 

Мы видели, что суммы Бореля приводят к инъективным гомоморфиз- 
мам дифференциальных алгебр и поэтому превосходно согласуются с вы- 
числениями с помощью расходящихся рядов. То же самое можно показать 
для суммы Стильтьеса (Ж.Экалль недавно замечательно обобщил соответ- 
ствующую идею под названием усредненное суммирование Бореля-Лапласа 
(см. также [Di], гл. 1, E., стр. 8 и [МВЗ], стр. 358). 

Этот результат тем более удивителен, что он не очевиден с перво- 
го взгляда. Действительно, рассмотрим, что происходит при разложении 
функции \/1 + т в сходящийся ряд в начале координат: 


1 Ly, Ш Oo 4 
Lea ee a ee a eee 
2 8 16 128 
Если мы хотим суммировать этот ряд при т = —2, то есть 
т. 
ПО 9: 
7. 2 8 
то одним из естественных способов было бы использование аналитического 
продолжения, избегая сингулярности при т = —1. Но тогда получатся два 


Половина разрывности в форме возникает, когда мы достигаем луча Стокса, а половина — 
когда покидаем его с другой стороны. 
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естественных решения, соответствующие двум ветвям функции \/1- 2; 
получаем +2. Ясно, что здесь нет естественного действительного решения. 
(Если бы f(a) была естественной действительной суммой для т < —1, мы 
бы получили f(x)? =1+х < 0, что невозможно.) 

В первом случае две суммы Бореля также соответствуют в некотором 
смысле (к которому мы вернемся в разделе 3.5) двум ветвям одной функ- 
ции (эта идея появляется уже у Стокса: аналогично смене знака в радика- 
ле |Sto2]), но смена ветви (явление Стокса) — в этом случае — происходит 
в параметрической группе. Действительно, автоморфизм дифференциаль- 
ной алгебры, соответствующий явлению Стокса, — это экспонента произ- 
водной, которая коммутирует с производной алгебры: посторонняя точечная 
производная Ж. Экалля. Подобная ситуация уже не имеет места для обыч- 
ной (алгебраической) смены ветви во втором случае (раздел 3.5). 


1.7. Пуанкаре и асимптотическая теория 


Классическая асимптотическая теория создана Пуанкаре. Последний 
разработал ее с намерением применить к аналитическим дифференциаль- 
ным уравнениям: основным его мотивом было придать смысл формальному 
решению такого дифференциального уравнения, записанному в виде рас- 
ходящегося ряда, то есть «воплотить» формальное решение в настоящее. 
Следует заметить, что принадлежащее Пуанкаре определение, которое мы 
сейчас приведем, не имеет ничего общего с его собственными замечаниями 
о свойствах изначально сходящихся, но затем расходящихся асимптотиче- 
ских рядов (то, что он называет сходимостью в смысле астрономов), как 
было бы естественно ожидать (и что привело к суммированию до наимень- 
шего члена в разделе 1.3). Сегодня известно, что в действительности это 
общее свойство формальных решений аналитических дифференциальных 
уравнений. 


Определение. Пусть У — открытый сектор с вершиной 0 на ком- 
плексной плоскости (или на римановой поверхности логарифма). Пусть f — 


oe OU) 

голоморфная Ha У функция и f = > ana” © C|[a]]. Говорят, что f асим- 
n=0 

птотически приближается рядом f на У, если для любого компактного 

подсектора W из У U {0} и любого целого п Е М существует действитель- 


ное Л/у/ в > 0 такое, что 
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для всех x € W. 


При выполнении этих условий мы обозначаем f ~ f или f= J(f). 

Множество голоморфных Ha У функций, допускающих асимптотиче- 
ское разложение в начале координат (снабженное операциями +, -, х?а/ат), 
является дифференциальной С-алгеброй, обозначаемой А(У). Имеем точ- 
ную последовательность дифференциальных алгебр 


0 A<°(V) > A(V) >, <] > 0. 


Сюръективность отображения J составляет теорему Бореля-Ритта (Воге]- 
Ritt): любой формальный ряд f может порождать” функцию f, голоморф- 
ную на У. К сожалению, это воплощение не единственно, функцию } можно 
выбирать (по модулю «пространства ошибок» AX<°(V), состоящего из го- 
ломорфных на У функций, бесконечно плоских в начале координат). Иначе 
говоря, у нас нет естественного вложения = f такого, что J(f) = f HeT 
теории суммирования в определенном выше смысле. Это один из наиболь- 
ших недостатков теории Пуанкаре (центральный недостаток классифика- 
ции Пуанкаре |Di| — отсутствие единственности для функции, представ- 
ленной асимптотическим разложением, составляет различие с суммами 
сходящихся рядов [О)). 

Применяя теорию Пуанкаре к аналитическим дифференциальным 
уравнениям, следует начать с формального ряда > как решения аналити- 


ческого уравнения 
(a, У, ..., y™) = 0, 


где G(x, %,..., Yn) — аналитическая функция п + 2 переменных 
и G(x, f,...,f™) = 0 и, воплотить Ё в истинное решение } 
(G(x, f,..., 1”) = 0). Тогда совокупность «возможных ошибок» (неуве- 


ренности относительно /) заметно уменьшается (для линейного уравнения 
она, очевидно, конечномерна, тогда как A<°(V) бесконечномерно), хотя 
в общем случае она, к сожалению, нетривиальна. За недостатком единствен- 
ности, мы, во всяком случае, имеем весьма удовлетворительный результат 
существования. Это фундаментальная теорема об асимптотических раз- 
ложениях. 


12 псагпег (франц.) — можно также перевести породить, материализовать. — Прим. ред. 
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Теорема. Пусть G(x, У, ..., Yn) — аналитическая функция п + 2 пе- 


ременных и ряд = С] || — формальное решение дифференциального урав- 
нения 


GO ee) SO (то есть G(x, a oe f™) ==); (1) 


Тогда существует действительное К > 0 такое, что для любого открыто- 
го сектора V с вершиной в начале координат, с углом раствора < п/К идо- 
статочно малым радиусом существует функция f, являющаяся решением 
дифференциального уравнения (1) и имеющая f своим асимптотическим 
разложением на V. 


В таком виде (то есть без каких бы то ни было ограничивающих ги- 
потез) этот результат был получен достаточно недавно, он принадлежит 
Рамису и Сибуя (Ramis et Sibuya) [К $1]. Первые частные случаи, установ- 
ленные Пуанкаре и [К 51], успешно улучшались многими авторами. 


2. Асимптотические разложения и суммируемость 


Сейчас мы вкратце обозначим, как можно преодолеть недостатки асим- 
птотической теории Пуанкаре, чтобы перейти к «точной асимптотической 
теории». Этапы таковы: 

® асимптотические разложения Жевре, 


e К-суммируемость, 


® мультисуммируемость 3. 


2.1. Асимптотические разложения Жевре 


Для начала сделаем несколько замечаний о явлениях, которые систе- 
матически встречаются при практическом использовании асимптотических 
разложений, но не приняты в расчет в теории Пуанкаре. 


(а) Можно констатировать, что большая часть встречающихся формаль- 
+00 
ных рядов >) апт” явно удовлетворяют условию вида 
n=0 


1 
lan| < C(n!)* A”, 


'S Многочисленные направления современных исследований развиваются в том же духе: 
теории рекурсий и ускорения сходимости Экалля [ЕП], [E2], [E3], работы русской школы: 
Ильяшенко, английской школы: гиперасимптотики [ВН]. 
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где С’, А, А > 0 — подходящие действительные числа. Для примера можно 
открыть одну из «библий», посвященных специальным функциям ([AS], 
[MOS], [Lu]), и подтвердить это наблюдение экспериментально. 

(6) Можно констатировать, что на практике неопределенность, вноси- 
мая в искомую функцию } с известным асимптотическим разложением", — 
это не только голоморфная функция, бесконечно гладкая на некотором сек- 
торе, но, точнее, голоморфная функция с экспоненциальным убыванием 


некоторого порядка k > 0: 
—а 


Race ee ei ет" (С. а>0 — подходящие действительные числа). 


(с) При изменении x номер наименьшего члена в общем случае бли- 


30K к номеру № = целой части + при подходящем выборе 6, К > 0 (для 
46; 


ряда Эйлера, например, А = 1 u 6 = 1). Это приводит нас к определению 


оо 
квазисуммирования до наименьшего члена: в качестве суммы ряда D> апт" 
n=0 
~ b 
берется число >) nx”, где N — целая часть —, параметры b и К выбра- 
n—0 LX 


ны заранее. Ha практике можно наблюдать, что при правильном выборе 
констант этот метод очень эффективен численно. 

Существует простое видоизменение асимптотической теории Пуанка- 
ре, превосходно объясняющее эти три замечания: это асимптотическая те- 
opus Жевре”?. 


Определение. Пусть А > 0 — действительное число. Пусть У — от- 
крытый сектор с вершиной 0 на комплексной плоскости (или на римановой 


поверхности логарифма). Пусть f — голоморфная на У функция, а } = 


+00 
= So anx” Е С |]. Говорят, что f асимптотически приближается no 
n=0 


Жевре порядка К рядом f на V, если для любого компактного подсекто- 
ра И’ из У U {0} и любого целого п € М существуют действительные 


14 Неопределенность, которая может возникнуть из-за погрешностей или из-за того, что есть 
несколько «естественных» решений задачи (двойственность). 

15Можно показать эквивалентность трех свойств (а), (b) и (с), если сформулировать их 
подходящим образом. Я отметил эквивалентность (a) < (5) в 1978 году; на эквивалентность 
(а) <= (с) мне несколько позже указал Б. Мальгранж. 
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Cw > Ои Aw > 0 такие, что 


п—1 


1 
ja|- |1 (=) — > сара] < C(nl) А", 


p=0 


для всех хе W. 


Множество голоморфных на У функций, допускающих асимптотиче- 


ское разложение по Жевре порядка s = 1 > ов начале координат (снаб- 


k; 
женное операциями +, -, х?а/4т), является дифференциальной С-алгеброй, 
обозначаемой A, (У). 
k 
Если У — открытый сектор с размыканием < л/К (малый сектор), то 


получаем точную последовательность дифференциальных алгебр 


0— ASV) = А: (У) > С] > 0. 
k k 
Сюръективность отображения J следует из теоремы Бореля-Ритта-Жевре 
(которую можно доказать с помощью «неполного преобразования Лапла- 
са»). 
Аналогом фундаментальной теоремы об асимптотических разложениях 


является фундаментальная теорема об ACUMNMOMUYECKUX разложениях NO 
Жевре [В $1]. 


Теорема. Пусть G(x, У, ..., Yn) — аналитическая функция n + 2 пе- 


ременных, а fe С]: — Формальный ряд Жевре порядка 7 являющийся 
k 


решением дифференциального уравнения 
G(x, y,..., 99) =0 (то есть G(x, f, ..., 9) =0). (1) 


Тогда существует действительное k' > 0 такое, что для любого открыто- 
го сектора У с вершиной в начале координат, размыканием < inf р 17) 
достаточно малого радиуса существует функция f — решение Oudde- 


ренциального уравнения (1), acumnmomuyecku приближаемая рядом f на V 


в смысле Жевре порядка i 


К 


Вся польза этой теоремы обнаруживается, если принять во внимание 
следующий результат. 
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Теорема (Мелле (Maillet) [M]). Пусть G(x, У, ..., Yn) — аналити- 
ческая функция п + 2 переменных, а f € С]|| — формальный pad, явля- 
ющийся решением дифференциального уравнения 


G(x, y,..., у) =0 (то есть G(x, f,..., 9) =0). (1) 


~_ 


Тогда существует действительное К > 0 такое, что формальный ряд f 
имеет тип Жевре порядка 1 /К. 


Мелле доказал эту теорему с помощью непосредственных оценок. В об- 
щем случае результат далек от оптимального (действительное число К не 
наименьшее из возможных). Мы вернемся к этому вопросу в разделе 3.1. 

Для подготовки к дальнейшему изложению (и чтобы объяснить взаи- 
мосвязь вышеприведенных замечаний (а) и (Ъ)) я собираюсь ввести «более 
геометрическую» концепцию ряда Жевре, используя понятие А-точной ква- 
зифункции. (На самом деле в литературе, на которую я ссылаюсь, в неко- 
торой степени используются когомологии пучков, но сейчас мне нужно 
передать основные идеи, оставаясь на элементарном техническом уровне.) 

Рассмотрим открытый сектор У на С (возможно, круг D* без од- 
ной точки) или на римановой поверхности логарифма. Используем по- 
крытия {У;}+сг на У. Предположим, что все эти покрытия конечны, pa- 
диус каждого У; совпадает с радиусом У, все тройные пересечения пусты. 
Предположим также, что Г = 1, т| и занумеруем V; по часовой стрелке. 
Обозначим Vi i414 = И; П Иди (если У = D*, то индексы считаются по 
модулю 77 + 1). Будем рассматривать голоморфные (-коцепи, совпадающие 
с последовательностями { f;}, где f; голоморфна на У; и соответствующие 
1-кограницы (1-коцепи), то есть последовательности { й, }, где №; = 1-Х. 


Определение. Пусть k > 0. Тогда АК-точная квазифункция на секто- 
ре У задается некоторой голоморфной 0-коцепью { f;}, ассоциированной 
с покрытием {У; ст, причем h; = fi41 — № убывают на У; ‚+1 экспонен- 
циально порядка К. Два таких набора ({}}; {У ет) и (195; ЗИ’; ел) 
отождествляются, если для каждого непустого пересечения У; ПИ’; раз- 
ность (f; — gj) убывает на этом пересечении экспоненциально порядка К. 


Иначе говоря, понятие К-точной квазифункции — это формализация 
понятия голоморфной функции, «известной с экспоненциальной точностью 
почти до порядка А». 
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Произвольная К-точная квазифункция называется ограниченной, если 
все f; ограничены. 

Покрывая проколотый диск D* секторами с размыканием < л/К 
и используя теорему Бореля-Ритта-Жевре, можно каждому ряду Жевре 
Fecal /к сопоставить единственную (по модулю проведенного отожде- 
ствления) К-точную квазифункцию, то есть его квазисумму (аналогичную 
функции, равной сумме сходящегося ряда). Если {f;} представляют ква- 
зисумму ряда f, тогда fj ~ f на У; в смысле Жевре порядка 1/k. Также 
как голоморфная функция, ограниченная на проколотом круге D*, является 
суммой некоторого сходящегося ряда (принцип несуществующих особен- 
ностей Римана), так и А-точная квазифункция, ограниченная на проколотом 
круге D*, является квазисуммой некоторого ряда Жевре порядка 1 /К. 

Это существенный результат, он доказывается как обычно, с помощью 
интеграла Коши ([Ка1], [Sil], [$12], [$13, th. 3.2]). Благодаря ему у нас есть 
«когомологический» метод доказательства того, что некоторый ряд имеет 
тип Жевре: это проверяется до «бесконечно гладкой поправки», чье экс- 
поненциальное убывание надо изучить. Как оказалось, этот метод мощнее 
обычных (то есть прямой оценки коэффициентов, теорем о неявных функ- 
циях и т. д.). Среди прочего, он допускает линеаризацию задачи (смотри 
часть 3). 

Основания асимптотической теории Жевре и теории К-суммируемости 
(раздел 2.2), найденные и развитые мной начиная с 1978 года, в действи- 
тельности датируются началом этого века и принадлежат английскому ма- 
тематику Г. Ватсону [ Wat2], [Wat3]. По-видимому, работы Ватсона не имели 
большого успеха, и его идеи были забыты. Еще сегодня можно найти от- 
голоски суровой критики, которой должен был подвергаться Ватсон, в двух 
наиболее известных (и лучших) книгах по асимптотической теории [Di], 
[О]: Дингл, упоминая об асимптотических разложениях Жевре и о введен- 
ной Ватсоном ассоциированной суммируемости, писал: 


Ценой заметного усложнения основная неточность определения Пу- 
анкаре может быть устранена... Было сказано достаточно для демон- 
страции всей сложности природы Этого (получившегося) определения... 
Доказательство среди этих рядов полных асимптотических разложений 
требует все более сильных и усиливающихся знаний (предположений). . ., 


16Однако следует отметить, что эти работы вызвали появление теории классов квазиана- 
литических функций Данжуа - Карлемана [Са] (также [Ne]), вновь возникшую в последних 
исследованиях Ж. Экалля о суммируемости расходящихся рядов: связанные функции. 
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чтобы реализовать обманчиво естественную идею, как это оказалось вна- 
yale, построить подходящий базис определений. Исключением являются 
простейшие асимптотические разложения, вычисляемые элементарными 
методами... 


Еще суровее — Олвер [O, стр. 543]: 


К сожалению, удовлетворительное определение полной пригодности 
недостижимо. Другой недостаток теории Ватсона — необходимость неко- 
торых свойств у остаточного члена, которые, по-видимому, достижимы, 
если только известна реалистичная граница остаточного члена. Так что 
эта теория, в основном, бесполезна. 


Мы надеемся, что следующая часть этих заметок убедит читателя в со- 
вершенной неточности этих суждений. (Очень удивляет ошибочная оценка 
у Дингла, который превосходно понимал важность экспоненциально малых 
ошибок в асимптотической теории и был так близок к правильной точке 
зрения!) 


2.2. К-суммируемость 


Пусть К > 0. Если У — открытый сектор с углом раствора < л/К 
(малый сектор), то, как мы видели, получается точная последовательность 
дифференциальных алгебр: 


0— AS*(V) + А. ж(У) > СШ 0. 


Совсем иная ситуация возникает, когда У — открытый сектор с раз- 
мыканием > л/К (большой сектор). Тогда получается следующая точная 
последовательность дифференциальных алгебр: 


0 Ayj(V) > С 0. 


В этом случае отображение J уже не сюръективно. Напротив, оно инъек- 
тивно, как показывает следующий результат. 


Теорема (Ватсон [Wat2]). Густь k > 0 действительно, а У — откры- 
тый сектор с вершиной в начале координат и раствором > п/Ё. Пусть f — 
голоморфная на У функция, убывающая на У экспоненциально порядка К. 
Тогда } тождественно равна нулю. 


2. АСИМПТОТИЧЕСКИЕ РАЗЛОЖЕНИЯ И СУММИРУЕМОСТЬ 4] 


Эта теорема — следствие теоремы Фрагмена-Линделефа (один их вари- 
антов принципа максимума). В случае большого сектора мы теряем суще- 
ствование, но получаем единственность суммы формального ряда Жевре. 
Это приводит к следующему понятию суммируемости. 


Определение. Пусть К > 0 действительно, а 4 — фиксированное на- 
правление. Формальный ряд f называется К-суммируемым по направле- 
нию 4, если существует функция f, голоморфная на некотором секторе V 
с биссектрисой 4 и размыканием > л/К, асимптотически приближаемая 
на У рядом f в смысле Жевре порядка 1 /k. 


Согласно теореме Ватсона при этих условиях функция } единственна 
(вблизи сектора определения У). Мы будем говорить, что f — это сум- 
ма ряда f по направлению d в смысле К-суммирования. Будем обозначать 


f = Spas, @ через Ст] : обозначим совокупность k-cyMMHpyeMBIX по 
р 
направлению 4 рядов. Сразу проверяем, что C{x}, — дифференциальная 
k ) 


подалгебра С]]т]] | ‚ а вложение 


2 |= 


инъективный гомоморфизм дифференциальных алгебр. (Здесь Ag обозна- 
чает алгебру ростков голоморфных функций на открытых секторах с бис- 
сектрисой 4, произвольным размыканием и радиусом.) 

Это определение К-суммируемости приемлемо (и очень полезно, на- 
пример, для проверки К-суммируемости ряда, являющегося формальным 
решением функционального аналитического уравнения), но оно не позво- 
ляет явно вычислять сумму. Для этого используется эквивалентное опреде- 
ление (эквивалентность легко доказать). 


Определение. Пусть k > 0 действительно, а d — фиксированое направ- 
ление. Формальный ряд f называется К-суммируемым no направлению d, 
если его формальное преобразование Бореля ¢ = В»! уровня Ё сходится 
и его (обычная) сумма аналитически продолжается голоморфной функцией 
(уже обозначенной) ф и возрастает экспоненциально порядка более чем А 
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в открытом секторе с биссектрисой 4. При этих условиях говорят, что 


в_- 2 
f(z) = в [96% Е а a) 
d 


это сумма ряда f no направлению d в смысле А-суммируемости. 


Оператор В»: С] — С]=]] был определен в разделе 1.5. Отметим, 
что оператор Лапласа Сь.а уровня К, введенный этим определением, задан 
на большей области, чем одноименный оператор, введенный в 1.5. Имеем 


Эна = Сва.абВь. 


Для К = 1 получаем метод суммирования Бореля (в действительности, 
изначальный метод Бореля имеет несколько более общий вид). 

Также К-суммируемый ряд можно суммировать по направлению 4, ис- 
пользуя один из абелевых методов, суммирование Харди-Жюрка (без явно- 
го использования параметра К). 


Теорема (Жюрка [Л]). Пусть К > 0, а а — исходящее из начала ко- 


ординат направление. Пусть f Е С] = У’. anx” — формальный ряд, 
п=0 

К-суммируемый по направлению d. Его сумма f по направлению d ана- 

литически продолжается вдоль максимального открытого интервала Эа, 

отложенного вдоль 4, функцией, которая тоже обозначается f. Если хо — 


точка U3 Yq, то предел 


+00 
lim (ло tay + aor? ae ) о 
t— 
t>0 n=3 


существует и равен f (xo). 


Легко видеть, что наша формулировка точнее, чем в [J]. Этот результат 
можно доказать — используя результаты автора этих заметок, — обратившись 
к случаю, когда коэффициенты a, А-суммируемого ряда »`аилх” являются 
моментами некоторой функции ф, определенной на [0, а] C В" и убыва- 
ющей в начале координат экспоненциально!” порядка k: 


"Eon К = 1, а функция ф принимает положительные действительные значения, сумму 
Бореля можно вычислить, преобразовав степенной ряд »`аьх?” в непрерывную сходящуюся 
дробь (Лагерр, Стильтьес, Борель [Bol], гл. 2, Падэ). 
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On = ] ФЕ Па. 
0 


Готовясь к изучению мультисуммируемости, мы переформулируем по- 
нятие К-суммируемого ряда немного более «геометрическим» способом, 
с помощью понятия квазифункции. 


Определение. Пусть А > 5 ufeé С] 1/к. Ассоциированная с f k-no- 


следовательность no направлению d — это последовательность (f1, }2), 
где fo — квазифункция, К-точная квазисумма ряда i. a fi — функция, го- 
ломорфная на секторе У на С с биссектрисой d и размыканием > л/К; 
функция f; как К-точная квазифункция равна ограничению fo на У. 


По теореме Ватсона, такая последовательность единственна, если она 
существует. Сразу проверяем, что такая К-последовательность существует 
тогда и только тогда, когда ряд f К-суммируем по направлению d; в этом 
случае }1 равна сумме ряда по направлению d B смысле К-суммируемости. 


2.3. Мультисуммируемость 


После того как Эмиль Борель открыл борелевскую суммируемость, 
а Леруа [Le] и Неванлинна [Ne] вскоре обобщили ее в понятие А-сумми- 
руемости, многие математики пытались доказать, что формальные ряды — 
решения дифференциальных алгебраических уравнений — всегда К-сумми- 
руемы. Уже в 1903 году эту задачу изучал Мелле [М]. К сожалению, как 
я показал в 1979 году в [Ва2 |, это неверно, причем по причине, в сущности, 
достаточно очевидной, так что удивительно, почему подобное замечание не 
было сделано раньше. 

Интуитивно понятно, что когда К = К', методы К-суммирования 
и К'-суммирования «достаточно далеки друг от друга» и почти несравнимы: 
если k < k’, то асимптотические оценки, нужные для ('-суммируемости, бу- 
дут строже, зато необходимое размыкание сектора (больше, чем л/К’) — 
меньше (требования к размыканию сектора менее строгие). Интуиция под- 
тверждается следующей «тауберовой» теоремой. 


Теорема [Ва7|. Пусть К и К' — два действительных числа, причем 
К’ > k > 0. Тогда 
Cllel], пс] 1 = Cla}. 
k к’ 
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Теперь, желая построить искомый контрпример, мы приходим к рас- 
смотрению суммы { формального К1-суммируемого решения ia линейного 
алгебраического дифференциального оператора 01 и формального А2-сум- 
мируемого решения ie линейного алгебраического дифференциального опе- 
ратора 15. Эта сумма будет решением линейного алгебраического диффе- 
ренциального оператора D. 

Наиболее простой пример получается, если взять в качестве fi ряд 
Эйлера, а в качестве fo — ряд Эйлера, в котором x заменен Ha 17. Тогда ft 
будет 1-суммируемый, В — 2-суммируемый ряд, но сумма 


f=At+h, К =У Сия + У (пума 
n=0 n=0 


не будет А-суммируема ни для одного действительного К > 0. Более того, 
сумма f — решение следующих линейных алгебраических дифференциаль- 
ных уравнений [RS], 3.75, 7.6]: 


e Dy = 0; 


о. Ш eee ee. a 2 | 
О = {2 (2 т) 7 2 x (22° — 5x 4) + 2(a x +2)}; 


e 2°(2—x)y” — 22(213 — 512 — 4)y’ + 2(12 — 2 + 2)y = - 321 + 1023 + 


Если нам нужен метод, пригодный для суммирования формальных pe- 
шений дифференциальных алгебраических уравнений (даже только линей- 
ных), то, как можно видеть, кажется разумным искать метод, позволяющий 
суммировать конечные суммы К-суммируемых рядов (по одному и тому 
же фиксированному направлению d) при разных К. Очевидно, можно поду- 
мать, что для f= a + fo сумму ряда f можно определить как сумму сумм 
рядов fi И р Но тут же возникает затруднение: надо доказать, что получен- 
ная таким образом сумма f не зависит от разложения (то есть от выбора fi 
И fa). Более того, мы желаем работать с алгеброй, поэтому нам надо изу- 
чить, как раскладывается в сумму произведение. Наконец, если мы хотим 
применить все это к дифференциальным уравнениям (линейным или нет), 
то нам нужно будет выяснить, как раскладывается в сумму А-суммируемых 
рядов (по одному и тому же направлению) формальный ряд, являющийся 
решением такого уравнения. 
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По прошествии многих лет — и благодаря усилиям многих математи- 
KOB — Ha все эти вопросы были даны удовлетворительные ответы. Проще 
всего определить мультисуммируемость, используя следующий результат. 


Теорема. //усть 4 — исходящее из начала координат направление. 
Сумма 


комплексных векторных подпространств из С |т|| является дифферен- 
циальной подалгеброй алгебры (С[[з||, х?4а/4х). Более того, существует 
единственный инъективный гомоморфизм дифференциальных алгебр 


oe C{t}1 > Ad 


такой, что Ограничение Sa на каждое пространство С{т} 1 (при к > 5) 
К ) 
совпадает с оператором К-суммирования по направлению d: 
Skid: С{т} 1... Ад. 
ard 
k 
Заменяя в формулировке теоремы переменную XY на одну из ветвей 
1 


переменной x™ (где т € №), получаем инъективный гомоморфизм диф- 
ференциальных алгебр (тоже обозначаемый 5х.а): 


Обозначим 
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которая содержит алгебры С{т} 1 „ для всех К > 0. Действительно, имеем 
т; 
равенство 
i 
Сет}: ПС] = С} a 
К’ mk’ 


Получаем инъективный гомоморфизм дифференциальных алгебр 
Oma Cix}..4 —> Ад, 


ограничение которого на каждую из алгебр C{x}, при К > 0 совпадает 
с; 
с оператором К-суммирования по направлению 4. 
Формальный ряд i принадлежащий C{x}..q, называется мультисум- 
мируемым по направлению 4, если для некоторого подходящего т Е N 


имеем ‚ 
i 
ТЕ >. С{хт т 
t= a 


Определение. Пусть ki, ..., К, — действительные числа такие, что 
ky > А >... > К, > 0, ad — исходящее из начала координат направление. 
Если f принадлежит 


r 1 
: | Cla™ } т 
i=1 ki’ 
то говорят, что ряд f (К1,...,К,)-суммируемый по направлению 4, и обо- 
значают f € С{т} 1 i <3 
5 a eee 
ky к, 


Существуют два достаточно различных доказательства только что 
сформулированной теоремы. Первое использует теорию ускорения Ж. Экал- 
ля [E4] и обобщает метод приближения интегральной формулой и ана- 
литическое продолжение (Бореля-Лапласа) для К-суммируемости. Это до- 
казательство подробно изложено в [MR3]. Второе, более геометрическое, 
принадлежит Б. Мальгранжу и Ж.—П.Рамису [МаК], оно существенно ис- 
пользует понятие экспоненциально малой поправки. Грубо говоря, метод 
таков. 

Пусть Аз, А2,..., К, — действительные числа такие, что ky > ko > 

. > К, > 0, a d — исходящее из начала координат направление. Пусть 
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ЕС тийИс--- СУ, — вложенные открытые секторы с одинаковым 


ky 

радиусом и биссектрисой. Предположим, что угол раствора V; > п/К.. 

Ассоциированная с } (ki, ..., К, )-последовательность задается по- 
следовательностью (}1,..., fr, fr+i), ге fi — функция, голоморфная 
на секторе Vj, f,11 — К,-точная квазифункция, определенная рядом f, 
а fo,..., fry — соответственно А1,..., К,_1-точные квазифункции, опре- 
деленные на секторах Vo, ..., И, соответственно; ограничение +1 на V; 
совпадает с А;-точной квазифункцией, ассоциированной с А;_1-точной ква- 
зифункцией f; для =1,..., г (по определению считаем Ко = +00). 


Относительная квазианалитическая теорема, которую мы сформулиру- 
ем дальше, позволяет увидеть, что если ассоциированная с f (ki, ..., Ky )-T0- 
следовательность существует, то она единственна. При этих условиях гово- 
рят, что f мультисуммируем по направлению d, a |1 — его сумма по направ- 
лению 4. Тогда легко проверить, что нь 1 — инъективный гомоморфизм 
дифференциальных алгебр ({1 допускает f в качестве асимптотического 
разложения в начале координат). После некоторых усилий получаем опи- 
санные выше свойства разложений в сумму А-суммируемых рядов. 


Теорема (Б. Мальгранж [Ма1]). Пусть К’ > k > 0. Пусть У — om- 
крытый сектор на комплексной плоскости (или на римановой поверхности 
логарифма) с вершиной в начале координат и раствором > п/ЁК. Пусть f — 
k-mounaa квазифункция на У, убывающая на У экспоненциально порядка К' 
(то есть } представима некоторой 0-коцепью { |+}, где р, убывают на V; 
экспоненциально порядка К, а +1 убывает на У; +1 экспоненциально по- 
рядка К). Тогда f — нулевая К-точная квазифункция на У (то есть все f; 
убывают на У; экспоненциально порядка К). 


Эта теорема обобщает сформулированную в 2.2 теорему Ватсона с ана- 
логичной формулировкой, в которой квазифункции заменены на функции. 
Относительная квазианалитическая теорема Мальгранжа эквивалентна «та- 
уберовой» теореме Ж. Мартине и Ж.-П. Рамиса [[MR2], chap. 2, prop. 4.3]. 


kk’ 
(k—k’) 
координат направление. Если формальный ряд fe Cl[z]]1/~ k’-cymuupyem 


по направлению d, mo он К-суммируем no любому направлению сектора 
с биссектрисой 4 и раствором т /К. 


Теорема. ГГусть 0 < hk’ < Е k= а 4 — исходящее из начала 
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Очень близкое к этой теореме утверждение получили еще Г.Х. Хар- 
ди [НЗ] (при т. = 2) и его ученик Гуд (Good) [Goo]. 


По поводу природы тауберовых теорем Г.Х. Харди писал в [Hl], 6.1, 
р. 121]: 

Существует другой, менее очевидный предел для эффективности этих 
методов и всего, что было доказано полезного. Ни одним методом нельзя 
суммировать ряд, расходящийся слишком быстро; но также нельзя сум- 
мировать расходящийся ряд со слишком медленной расходимостью". Te- 
оремы, воплощающие этот принцип, принадлежат... к так называемому 
тауберову классу. В них утверждается, что если ряд суммируем (Р) и удо- 
влетворяет еще какому-то условию Кр (которое меняется в зависимости 
om Р, но, во всяком случае, влечет за собой некоторую замедленность 
возможной расходимости), то он сходится. 


В некотором смысле этот принцип был обобщен позже, согласно идеям 
Харди и Литлвуда. Процитируем Гуда [Goo, р. 145], который сравнивает два 
метода суммируемости f и g: 


Итак, очень часто мы получаем принцип суммируемости Харди-Литл- 
вуда: если метод } имеет более широкую область приложения, чем ц, 
и если g приложим, а f эффективен, тогда д эффективен. 


Сформулированная выше тауберова теорема хорошо иллюстрирует 
этот последний принцип. 

Недостаток приближения с помощью только что очерченного нами ме- 
тода мультисуммируемости состоит в том, что у нас нет явного удобного 
способа для вычисления суммы. Можно получить такой способ с помощью 
ускорения, которое мы здесь описывать не будем (можно обратиться к ввод- 
ной статье в [LR2]). Еще можно использовать метод итераций преобразо- 
ваний Лапласа (разных уровней), принадлежащий Бальсеру (Balser) [Ba2]. 
Он действует так. 


Пусть kj >... > К, > 0. Определим строго положительные действи- 
тельные числа А1,..., ky условиями kp’ = Кик = k>*— ko, при 
1=2,..., Г так, что 

1 1 1 
SS ree ATID: о 
ks ky ky p ’ ’ 


Обозначим через .а аналитическое продолжение по направлению 4. 


18Выделено Харди. 
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Теорема. //усть 4 — исходящее из начала координат направление, 
ky >... >. > Ouk,, ..., kp определены выше. Для f € С] следующие 
условия эквивалентны: 

(i) ряд a [В ‚.› Кг-суммируем по направлению 1); 


a 


(ii) ряд By... Bx, f сходится, и для всех 1 =1,..., 2 функция 
аСк..а... Ск абВь....Вь f 


голоморфна и возрастает экспоненциально порядка К;_1 на подходящем 
секторе с биссектрисой 4. . ee 

При выполнении этих условий функция Сь..а...Ск,..аоВБе....Вы} су- 
ществует и равна сумме (в смысле мультисуммируемости) ряда f no Ha- 
правлению 4. 


Этот метод допускает численные приложения [Th2]. Другой явный Me- 
тод суммирования мультисуммируемого ряда — элегантный, но менее под- 
ходящий для численных приложений — задается следующим результатом 
Жюрка [J]!’. 


Теорема. Пусть 4 — исходящее из начала координат направление 
ze 00 
uf € Cla]] = So anx” — формальный ряд, мультисуммируемый no на- 
=0 


правлению d. Его сумма f по направлению d аналитически продолжает- 
ся вдоль максимального открытого интервала Yq, отложенного вдоль 4, 
функцией, тоже обозначаемой f. 

Если хо — точка на Эа, то предел 


+00 
lim (о Е aor as | И. 
t— 
t>0 n=3 


существует и равен f (xo). 


Иначе говоря, мультисуммируемость включает суммируемость абеле- 
вым методом Харди-Жюрка. 

Заметим, что для разработки явных методов суммирования, осно- 
ванных на мультисуммируемости (суммируемость ускорением или итера- 
циями преобразований Лапласа), необходимо знать критические уровни 


1ЭВ действительности доказанный Жюрка результат менее точен. Нужно обобщить его до- 
казательство, используя [Мак |. 
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К1,..., kp. Напротив, для суммирования абелевым методом Харди-Жюр- 
ка их знать не надо: последний метод «мощнее» (он суммирует больше 
расходящихся рядов). Мы встретились с ситуацией, аналогичной суммиро- 
ванию сходящегося ряда в его звезде Миттаг-Леффлера: когда абелев метод 
Линделефа позволяет находить сумму во всей звезде, метод Бореля-Лап- 
ласа (Бь, Ск) требует выбора параметра К в зависимости от точки звезды, 
в которой мы желаем найти сумму. За мощность метода надо платить чис- 
ленной неэффективностью. 

Такова иллюстрация общего феномена, открытого Харди и Литлвудом: 
чем мощнее метод суммирования, тем он чувствительнее (тоньше). Вновь 
процитируем Харди [НЗ], р. 153: 


Мы с Литлвудом часто настаивали на общем принципе, который 
трудно сформулировать точно, но который в общих чертах можно обо- 
значить так: чувствительность метода суммирования стремится быть 
обратно пропорциональной его мощности. 


cP Расходящиеся ряды и динамические системы 


3.1. Формальные решения дифференциальных уравнений 


Здесь я не буду говорить особо о линейном случае. По поводу 
численных аспектов можно обратиться к обзору М. Лодэ-Ришо (Loday— 
Richaud) [LR1] ик [Th1] u [Th2]. 

Сейчас я вернусь к теореме Мелле, сформулированной в 2.1. Целью 
Мелле было просуммировать формальные ряды, являющиеся решениями 
алгебраических дифференциальных уравнений. Он думал, что здесь можно 
добиться успеха с помощью К-суммируемости; установленное им свойство 
Жевре было бы необходимым условием. Мы видели, что в действительности 
этот план сам по себе нерационален, и оказалось необходимым прибегнуть 
к мультисуммируемости. Сейчас мы увидим, что ее достаточно: любой фор- 
мальный ряд, являющийся решением аналитического дифференциального 
уравнения (линейного или нет), мультисуммируем по всем направлениям, 
кроме разве что некоторого конечного числа. 

Для начала улучшим условия Жевре в теореме Мелле. Простейший 
случай — линейный, с которого мы и начнем. 


Определение. Будем говорить, что формальный ряд f имеет mun Жев- 
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ре точного порядка 1/k, если он имеет порядок 1/k, и не существует дей- 
ствительного k’ > К такого, чтобы ряд имел порядок 1 /К. 


С каждым ростком линейного аналитического дифференциального 
уравнения 
Dy = any™ +... Наву = 0, 


в начале координат свяжем многоугольник Ньютона N(D) [Ral]. Имеет 
место следующий результат. 


Теорема. Густь Dy = 0 — росток аналитического дифференциально- 
го уравнения в начале координат комплексной плоскости. Если f € С] — 
формальное решение этого уравнения, то f либо сходится, либо имеет 
mun Жевре точного порядка 1 /К, где К — одна из строго положительных 
сторон многоугольника Ньютона N(D) для О. 


Этот результат получен О. Перроном в случае линейного алгебраиче- 
ского дифференциального уравнения и мной — в аналитическом случае [Pel], 
[Ral, 5]. (В линейном аналитическом случае оценки Мелле были улучшены 
в работе Джингольда [G1].) 

Перейдем к нелинейному случаю. Рассмотрим аналитическое диффе- 
ренциальное уравнение 


G(x, y,...,y™) =0. (1) 


Предположим, что это уравнение допускает формальное решение f Е С [2%]. 
Мы можем сопоставить паре (С, Г) многоугольник Ньютона N(G, }) (мно- 


гоугольник Ньютона функции С вдоль 1). Можно получить точные оценки 
в теореме Мелле. 


Теорема. Пусть G(x, Yi, ..., Yn) — аналитическая функция n+ 2 
переменных и ЕС] |] — формальный ряд, являющийся решением диффе- 
ренциального уравнения 


О у”) ==:(), (1) 


Тогда } сходится или имеет тип Жевре точного порядка 1 / К, где К — одна 


— 


из строго положительных сторон многоугольника Ньютона N(G, f). 


Сначала Б. Мальгранж показал, что формальное решение всегда имеет 
тип Жевре порядка 1/k’, если k’ — наименьшая сторона N(G, }) [Ma3]. 
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Сформулированный выше результат принадлежит И. Сибуя [Sil], который 
доказал его с точки зрения когомологий, оценивая с помощью линеаризации 
1-коцепь экспоненциальных ошибок, ассоциированных с квазифункцией, 
определенной рядом р 

Предыдущие результаты улучшаются с помощью мультисуммируемо- 
CTH. 


Определение. Пусть А1 > ko >... > kp > 0 — действительные числа. 
Мы будем говорить, что РЕ С] является (k1, ..., К,)-суммируемой, если 
она (ki, ..., К,)-суммируема по всем направлениям, кроме, может быть, 
некоторого конечного числа. 


Теорема. Пусть G(x, У1,..., Yn) — аналитическая функция n+ 2 


переменных и f © С]т|] — формальный ряд, являющийся решением диффе- 
ренциального уравнения 


G(x, y,..., у) =0. (1) 
Пусть ky > kg >... > kp > 0 — строго положительные стороны много- 
угольника Ньютона N(D, 1). Тогда f (ki, ..., К,)-суммируема. 


В линейном случае первое доказательство этого результата принадле- 
жит автору этих заметок [Ка2], [MR3]; другие доказательства можно найти 
в [ВВК$] и [Mar]. В нелинейном случае этот результат недавно доказал 
Брааксма (Braaksma), использовавший подход Ж. Экалля. В настоящее Bpe- 
мя готовится доказательство Рамиса и Сибуя в стиле [В $1]. 


3.2. Нормальные формы дифференциальных уравнений 
и диффеоморфизмы. 


Сначала рассмотрим случай локальных голоморфных диффеоморфиз- 
мов в начале координат комплексной плоскости. Такой диффеоморфизм 


f: (С, 0) — (С, 0) 


определен сходящимся рядом 
оо 


f(z) = + У ane” Е xC{z}, 


n=1 


toe A € С”. 
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Пуанкаре доказал, что при |Л| ~ 1 такой диффеоморфизм аналитиче- 
ски сопряжен со своим касательным линейным отображением fo: х ++ Ах 
в начале координат. Иначе говоря, существует аналитическая замена пе- 
ременных x = y(t), удовлетворяющая условию w’(0) 0 и такая, что 


f= vofooy' или fop=po fo. 


e Если |Л| = 1, то говорят, что А лежит в области Пуанкаре, иначе 
говорят, что А лежит в области Зигеля. 


e Если Л — корень из единицы, то говорят, что это резонансный диффео- 
морфизм. 


® Если |Л| = 1, но \ — не корень из единицы, то можно доказать, что } ли- 
неаризуема формально (существует формальная замена переменных u, 
сопрягающая f u fo), HO He всегда линеаризуема аналитически (ряд w 
не обязательно сходится). Мы больше He будем задерживаться Ha этом 
сложном случае, его изучение связано с задачами об аппроксимации 
действительных чисел рациональными. 


Нас интересует резонансный случай. Диффеоморфизм f в общем слу- 
чае уже не будет формально линеаризуемым, и первая из возникающих 
здесь задач состоит в формальной классификации резонансных аналитиче- 
ских диффеоморфизмов. Чтобы упростить изложение, отныне мы ограни- 
чимся диффеоморфизмами f, касательными к тождественному: 


со 
f(z) =т+ У. тт" € Ст}. 
1 


Каждый класс формальной эквивалентности (то есть по модулю фор- 
мальной замены переменных) представляется (сходящейся) нормальной 
формой, выбранной произвольно (ищут «самую простую из возмож- 
ных форм»). В качестве нормальных форм удобно выбирать множество 
Тв.к.х = ехр(Хв.к.л) с параметрами (3, К, Л) Е C x №" х С, где Хвьл — 


векторное поле 


ak 


1+ xk 


Можно легко свести все к случаю В = 2i7 и обозначить Хожк,х = Хкл. 
Чтобы еще больше упростить изложение, мы ограничимся случаем, когда 


Хв,к,л = В 
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нормальная форма равна 


fio(x) = —= 


1 — 2172 


с X10 = 2inx7d/dzx. Гомографическая замена координат T= т на римано- 


вой сфере сопрягает нормальную форму }1.0 с переносом Т(=) = z — 2i7. 
Заметим, что функция е/. постоянна на орбитах fio. Она позволяет отож- 
дествить пространство орбит для f1.9 с римановой сферой с разрезом от 0 
до со, то есть с С* (которое топологически эквивалентно цилиндру 51 x R). 
В дальнейшем будем обозначать fy = f1.0. Имеем 


fo(x) =х-— пд? — 4243 +.... 
Любой локальный диффеоморфизм, голоморфный форме 
f(x) = 2 — п? — 4п? 13 + O(c'), 


формально сопряжен с fo (HO, в общем случае, не сопряжен аналитиче- 
ски). Точнее, существует единственный формальный | диффеоморфизм w, 
касательный к тождественному и такой, что фо = w о fo. Этот диффео- 
морфизм «в общем» расходится, но важно то, что он 1-суммируем по всем 
направлениям, кроме действительных полуосей R* и В `. Суммы фа диф- 
феоморфизма w по разным направлениям d суммируются, когда d меняется, 
не пересекая действительной оси. Итак, в конечном счете получаем две сум- 
мы wT uw (верхнюю и нижнюю) диффеоморфизма w. Они определены 
на соответствующих секторах с размыканием 37. Пересечение этих секто- 
ров состоит из двух секторов с размыканием 7, содержащихся в Rex > 0 
и Rex < 0 соответственно. В каждом из этих секторов сумму можно опре- 
делить двояко, здесь и возникает феномен Стокса. 

Обратимся теперь к классификации голоморфных локальных диффео- 
морфизмов, формально сопряженных с fo, по аналитическому модулю экви- 
валентности (два диффеоморфизма эквивалентны, если они сопряжены ана- 
литически). Удивляет тот феномен, что фактор-пространство очень велико: 
оно бесконечномерно. Оно существенно параметризуется парами (Yo, Фо), 
где Yo и Yoo — локальные диффеоморфизмы римановой сферы, касательные 
к тождественному преобразованию в 0 и в OO соответственно, но, вместе 
с тем, совершенно произвольные. 

Действительно, локальный голоморфный диффеоморфизм f, формаль- 
но сопряженный с fo, существенно классифицируется своим пространством 
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орбит. Это пространство легко описать с помощью феномена Стокса. Оно 
представляет собой «связную (неразделенную?) комплексную кривую», по- 
лучаемую следующим образом. 

Каждая сумма wT ww диффеоморфизма 1, сопрягающего f и fo, поз- 
воляет отождествить некоторое подмножество пространства орбит } с про- 
странством орбит fo, то есть с цилиндром C* = Р"(С) — 10, co}. Если 
сравнить эти два отождествления с помощью двух явлений Стокса по на- 
правлениям R* и В , то получим две копии римановой сферы P!(C), 
склеенные в 0 ив CO аналитическими диффеоморфизмами фо и Yoo, ка- 
сательными к тождественному. Можно доказать (теорема синтеза), что эти 
диффеоморфизмы выбираются произвольно. 

Предыдущие результаты получены Т.Кимурой, Ж.Экаллем [E2], 
Б. Мальгранжем [Ма4] и Ворониным. 

Классифицировав резонансные диффеоморфизмы, можно перейти 
к аналогичной задаче для ростков в (0, 0) — начале координат С? — анали- 
тических дифференциальных уравнений вида w = Рау- Фатх = 0, где P 
и @ принадлежат C{x, у}. 1-струя функции и в начале координат имеет 
вид Jtw = Ах ау + py dx. Если Л, ИЕ С*иЛ/и ЕС - В, то мы находимся 
в области Пуанкаре и w аналитически сопряжено с J'w. Иначе рассуждения 
усложняются. Как и для диффеоморфизмов, нас будут интересовать только 
резонансные случаи. Это случаи, когда А или и (но не оба вместе) рав- 
ны нулю (вырожденный случай) или когда А, / не равны нулю и A/p Е QU 
(невырожденный случай). Геометрически (ссылаясь на действительный слу- 
чай) в вырожденном случае имеем узловую горловину?°, в невырожденном — 
резонансную горловину. 


Узловая горловина описывается уравнением x*tldy + Aydx +... = 0, 
глеке №илЕС. 

Резонансная горловина описывается как рт Ау + qydx+... = 0, где 
р, аЕ №. 


Как и в случае резонансных диффеоморфизмов, формальные нормаль- 
ные формы параметризуются конечным числом параметров, а замены пере- 
менных, сводящие резонансное уравнение к его нормальной форме, требуют 
использования рядов типа Жевре. Можно доказать, что эти ряды А-сумми- 
руемы. При фиксированной формальной нормальной форме резонансные 
дифференциальные уравнения аналитически существенно классифициру- 


20 поеи4-со| (франц.) 
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ются пространствами их листов (т.е. слоения). Это пространство описы- 
вается с помощью феномена Стокса. 

Действительно, классификация уравнений связана с классификацией 
диффеоморфизмов посредством голономии. Узловая горловина всегда до- 
пускает гладкий аналитический лист в качестве решения в начале координат 
(сильное многообразие). Резонансная горловина всегда допускает два таких 
решения (трансверсальных). Конечно, требуется удалить начало координат: 
тогда эти решения существуют в окрестности начала координат в проко- 
лотых кругах. Можно изобразить комплексную трансверсаль к одному из 
этих проколотых кругов, затем — простой разрез с началом на трансверса- 
ли, к проколотому кругу. Используя этот разрез, мы получаем перестановку 
листов, которую рассматриваем как диффеоморфизм на трансверсали: это 
диффеоморфизм голономии. В случае резонансной горловины классифика- 
ция уравнений отождествляется с классификацией их голономий. В случае 
узловой горловины мы не получим всех голономий. В качестве примера 
возьмем формальную нормальную форму узловой горловины: 


wo = т?ау + ydx = 0. 


ВЕБ. НЕЕ 
(1— 2415). 

Аналитические узловые горловины, формально сопряженные с Wo = 0, 
существенно классифицируются парами локальных диффеоморфизмов 
(аналитических и касательных к тождественному) (Yo, Yoo) в 0 ив со Ha 
римановой сфере (где фо произвольно, но фо — перенос). 

Эти результаты получены Ж. Мартине и Ж.-П. Рамисом [MRI], [MR4], 
[MR5] и [Ма4]. Они играют фундаментальную роль при ответе на некоторые 
предположения Р. Тома о голоморфных листах [Моч] и в недавнем решении 
задачи Дюлака (конечность количества предельных циклов для алгебраиче- 
ского дифференциального уравнения Рау + Фатх = 0 на действительной 
плоскости)". 


Ее голономия fo(x) = 


3.3. Сингулярные возмущения, запаздывание бифуркации и утки 


Как мы видели, теория асимптотических разложений Жевре очень эф- 
фективна при изучении существенных особенностей дифференциальных 
уравнений. Речь идет об особенностях переменных. Эта теория точно так же 


21Задача Дюлака — это подзадача второй части шестнадцатой проблемы Гильберта. 
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используется при изучении некоторых задач, требующих введения особен- 
ностей в параметре: задач о сингулярных возмущениях. Вот существенный 
результат в данном направлении. 


Теорема (И. Сибуя). Пусть п Ее № — фиксированное целое число. 
Рассмотрим дифференциальное уравнение вида 


а = F(x, в, у) = f(a, =) + A(a, e)y + У. ЕР. (1) 


|p|>2 


где с Е №", Е тЕ С, уе С", f, Еи все fp принимают значения из С", 
а А принимает значения из End(n, С). Предположим, что Е — аналити- 
ческая по (т, =, y) в окрестности (0, 0, 0), }(0, 0) = 0, а матрица А(0, 0) 
обратима. При этих условиях: 

(i) уравнение (1) допускает единственное формальное решение вида 


где An — функции, голоморфные в одном и том же круге D плоскости 
переменных т. 
(ii) При уменьшении круга D функции An удовлетворяют оценкам Жев- 
ре порядка 1/0: 
[а || < Си! ® Ап, 


при подходящих С, А > 0 (здесь ||an|| = зархер ||@n(x)|l). 


И. Сибуя дал когомологическое доказательство этого результата в [S12]. 
В действительности он доказал более точный результат, построив 0-ко- 
цепь } решений уравнения (1), где f; голоморфны по (x, =) на произве- 
дении D и сектора по =, допускают асимптотическое разложение в ряд f 
типа Жевре порядка 1/o по = равномерно по x из О. М. Каналис- Дюран 
(Canalis— Durand) дал доказательство этой теоремы при п = 1 с помощью 
прямых оценок [CD2]. С другой стороны, Шеффке (Schaffke) недавно по- 
лучил новое доказательство, обратившись к уточненной форме теоремы 
Мелле [Sc]. 

Следуя идее Ж. Мартине [Мак |, мы применим вышеуказанный резуль- 
тат к задаче об запаздывании бифуркации. Здесь будет удобно перенестись 
в рамки нестандартного анализа"”. 


22 Дальнейшее легко можно прочитать на эвристическом уровне. В действительности все 
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Основная идея — принадлежащая Ж. Мартине [MaR] — состоит в том, 
чтобы смоделировать классическую среду численного анализа: вычисли- 
тельная машина (возможно, это математик, снабженный бумагой и каран- 
дашом) обеспечивает некоторую ограниченную точность при вычислении 
значений функции. Возможность обращения с большими числами ограниче- 
на; слишком маленькое число считается нулем, слишком большое рассмат- 
ривается как бесконечность (переполнение). Эту ситуацию можно смодели- 
ровать, задав раз и навсегда бесконечно малое действительное число =’ > 0: 
мы «не увидим» комплексное число @ порядка =’ (иначе говоря, такое, что 
a/c’ ограничено, то есть не бесконечно велико); с другой стороны, @ не 
будет «показано», если его модуль не бесконечно мал по сравнению с 1/=' 
(то есть если ас’ не бесконечно мало). 

Тогда можно переделать предыдущее описание квазифункций, заме- 
HAA экспоненциально малые поправки на поправки порядка =’. Получаем 


словарь для двух точек зрения на одну и ту же область (на плоскость — 


k 
a/€" (со стан- 


нестандартную — комплексной переменной =), где функция е` 
дартным или ограниченным а) имеет порядок =”. 

Явление запаздывания бифуркации происходит от феномена экспонен- 
циального сжатия-взрыва траекторий. Простейшая из ситуаций, в которой 


можно наблюдать этот феномен, следующая. Рассмотрим систему диффе- 
ренциальных уравнений 
В о (1) 
и —&, 


где = > 0 — бесконечно малое фиксированное действительное число, а LU 
и x — действительные переменные. Если ненулевой параметр и фиксирован, 
получаем дифференциальное уравнение 


LS ie: 


Точка x = O — устойчивая стационарная при [1 < 0 и неустойчивая при 
и > 0. Если изменять LL, то при js = 0 устойчивость меняется на неустойчи- 
вость. Теперь рассмотрим решение (1), определенное начальными услови- 
ями 1 = Xo и и = Lo < 0. Обычное интегрирование показывает, что реше- 
нием здесь будет в 


т = тое 2 


обсуждается очень строго. Использована точка зрения Нельсона в ANS [DR]; необходимая 
величина из ANS... бесконечно мала. 
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Оно бесконечно мало (экспоненциально по =) при и Е [Uo, —по|. Оно 
«убывает почти вертикально» от (40, Xo) до (ро, 0), продолжается «экс- 
поненциально почти вдоль» действительной оси вплоть до (—/40, 0), затем 
«поднимается почти вертикально». Имеем (почти. ..) «вход» в точке (40, 0) 
и «выход» в точке (—ро, 0). 

Теперь рассмотрим систему (размерности р + 1): 


4 5 = F(z, д 


dt 
dt 


где = > 0 — бесконечно малое фиксированное комплексное число (или 
«малый параметр», по выбору), х — переменая из С?, Л — комплексная 
переменная, Ё`— голоморфная функция со значениями в СР. Из системы (2) 
можно вывести дифференциальное уравнение 


dx 
— = F(a, Л). 3 
Е = F(x, №) (3) 
Полагая = = 0, получаем «медленную кривую» F(x, A) = 0. Предположим, 
что эта кривая Со — график аналитической функции © = со(Л). 


Теперь предположим следующее: 


(1) Матрица Fo, 0) обратима. 


(11) Медленная кривая Со трансверсальна полю 


F(z, Nx a ek. (27) 

Если медленная кривая инвариантна в этом поле, мы получаем ситу- 
ацию, очень близкую к вышеописанной. Это — экспоненциальное сжатие- 
взрыв. Случай, когда выполняется условие (11), сложнее. Мы покажем, что 
здесь всегда есть экспоненциальное сжатие-взрыв (в выбранной подходя- 
щим образом окрестности начала координат), «сводя все к предыдущему 
случаю», то есть используя экспоненциально точную квазикривую (ква- 
зи)инвариантную в поле (2’). 

Эта квазикривая является «графиком» некоторой квазифункции — (ква- 
зи)решения уравнения (3), полученного следующим образом. По сформу- 
лированной выше теореме Сибуя уравнение (3) допускает формальное ре- 
шение (формальное по =, аналитическое по A): 


g(A, =) = со (Л) + У. Е”. 


n21 
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Более того, все Cyn — ограниченные голоморфизмы на одном и том же кру- 
ге Г) (стандартного) радиуса т > 0 на плоскости переменных Л и удовле- 
творяют неравенствам Жевре первого порядка: 


|Сп|| < CntA”, 


при подходящих С’, А > 0 (здесь ||сп|| = supyep ||Cn(2)||). 

Квазисумма этого ряда (полученная, например, с помощью неполного 
преобразования Лапласа) |-точна, она будет квазирешением (то есть реше- 
нием до экспоненциальных погрешностей близкого порядка, типа е—@/е ‚ где 
а — ограниченное положительное число). В действительности можно найти 
представление { д; } этого квазирешения, где д; — точные решения (3). Отсю- 
да легко делаем вывод (например, используя аргументы «экспоненциальной 
выпуклости», см. [ВСОО]). 

Первые математические результаты в задаче об опаздывании к бифур- 
кации принадлежат Нейштадту. Его метод — вариант вышеописанного ква- 
зисуммирования до наименьшего члена: осуществляется некоторое (боль- 
шое) число N замен переменных (где N — целая часть 1/=). Еще существу- 
ет достаточно геометрический подход к задаче, предложенный Ж.-Л. Калло 
(J.-L. Callot) в 1991 году. 

Вместо того чтобы брать медленную производную векторного поля, 
как в только что изученной нами задаче, можно взять медленную произ- 
водную от отображения итерации. Пусть везде = > 0 — бесконечно малое 
действительное число. Расмотрим отображение 


Fe: В >В; 


(4) 

(tz, A) (f(a, A), А-а). 
Предположим, что существует аналитическая кривая Co(A), состоящая из 
фиксированных точек. Можно найти медленную производную отображения 
Фейгенбаума: 


которое дает нам пример хаоса. Наблюдаем запаздывание бифуркации при 
«удвоении периодов». 

Эту ситуацию анализировали А. Фрюшар [Fl] и С. Бэсенс (С. Вае- 
sens) [Bae]. Методы [F1] следуют тому же направлению, что и Нейштадт; 
методы [Bae] используют аргументы Ж. Мартине для случая полей. 
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Для отображения (4) инвариантной кривой будет график функции 
Aw U(A, =), где Е(И(А, =), А) = U(A +6, =). Если 0, F(A, со(Л)) = 1, 


то существует единственная формальная инвариантная кривая 


где с„ — ограниченные голоморфизмы на одном и том же круге. В общем 
случае этот ряд расходится, но имеет тип Жевре 1 [Вае]. Как и в случае 
полей, получаем квазиинвариантную квазикривую. 

Теперь мы увидим, что если в анализе Ж. Мартине отказаться от усло- 
вия обратимости (1), то могут произойти очень интересные явления: появ- 
ляются утки. Для этого мы обсудим уравнение Ван дер Поля (в котором 
феномен утки был открыт в [BCDD]) с точки зрения Жевре. Здесь мы сле- 
дуем (высказывая оригинальную точку зрения на феномен уток) недавней 
работе М. Каналиса- Дюрана [CD1], [CD2]. 

Рассмотрим сингулярно возмущенную форму уравнения Ван дер По- 
ля (1920) с действительным бесконечно малым = > 0: 


=% + (x? —- Пет =0. (5) 
Полагая и = =ЕХ + 328 — т, переходим к илоскости Льенара. Получаем 


и 
ELS = =a ee, 
3 (6) 


и=-х. 


Векторное поле ((u — 1/32° + x)/e, —x) «медленно-быстрое»: его горизон- 
тальная составляющая бесконечно велика за исключением медленной кри- 


вой — кубической параболы и = 52 — т, на которой она обращается в нуль. 
Часть кубической параболы между точками В (—1, =) и A(1, — =) — оттал- 


кивающая, остальная часть — притягивающая. 

Можно показать, что система (5) имеет предельный цикл. Этот цикл 
«медленно-быстрый». Он бесконечно близок к стандартному циклу, обра- 
зованному двумя дугами медленной кривой (заканчивающимися в Аи В) 
и двумя горизонтальными отрезками (исходящими из Аи В), которые на- 
зывают его тенью. Кроме того, у системы есть неустойчивая неподвижная 
точка (0,0). 
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Теперь рассмотрим уравнение Ван дер Поля со вторым членом 
= +(x? —1l)t+a=a, (7) 


где а — действительный параметр. Будем изменять этот параметр от 0 до 2 
и следить за ситуацией на плоскости Льенара. Неподвижная точка 5(а) пе- 
ремещается по медленной кубической параболе: 5(0) = (0, 0) и 5 (1) = А. 


e для 0 <а < 1 эта точка неустойчива; 


e для а > | она устойчива, но предельный цикл исчезает. 


Значение а = 1 — это бифуркация Хопда. Задача в TOM, чтобы понять, 
как может проявиться эта бифуркация, если принять во внимание медлен- 
но-быстрый характер поля. Наблюдая за изменением тени предельного цик- 
ла в момент его исчезновения, мы заключаем, что при некоторых частных 
значениях эта тень должна на момент вытянуться вдоль отталкивающей ча- 
сти медленной кривой (от Ак В). По определению, соответствующие значе- 
ния параметра называются значения с утками, а соответствующие циклы — 
утками (с головой или без головы, в зависимости от того, находится точка В 
в их тени или вне ее). 

Чтобы получить значение с уткой, очевидно, необходимо, чтобы а было 
меньше | и бесконечно близко к 1, но хотелось бы знать о нем больше. Тогда 
можно доказать, что все значения а с уткой имеют одно и то же разложение 


по с-тени: 
1+ > ОЕ 
n21 


Это приводит к мысли, что «жизнь уток коротка». Действительно, их жизнь 
экспоненциально коротка: если а1 И ag — значения с уткой, то |а1 — 
—ag|| < e~°/* для подходящего ограниченного 6. Экспериментально можно 
установить, что: 


e При= = x уток можно наблюдать примерно между 0,9934909 


и 0,99349 15. 
e При = = a уток можно наблюдать примерно между 0,9987404512 


и 0,9987404513, то есть на промежутке длиной 10719... Так что очень 
сложно «поймать» уток, даже если известно их разложение по =-тени. 
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Все это привело меня где-то в 1980 году к гипотезе, что разложение по 
=-тени значений с уткой для уравнения Ван дер Поля имеет тип Жевре |. 
Если оставить в стороне теоретические вопросы, положительный ответ на 
эту гипотезу был интересен мне для получения численного метода охоты 
на уток: если ряд имеет тип Жевре, то, суммируя его с помощью неполно- 
го преобразования Лапласа (или квазисуммирования до наименьшего чле- 
на), найдем значение, определенное до экспоненциально малой поправки 
порядка почти =. Итак, ошибка — величина того же порядка, что и про- 
должительность существования искомого значения с уткой. Эта гипотеза 
(с ее численным приложением) сначала была проверена экспериментально, 
а недавно ее доказал М. Каналис-Дюран [СО1]|, [CD2], используя среди 
прочего улучшение вычислений при рекуррентной зависимости коэффици- 
ентов с„, изложенное в [ZS]. 

М. Каналис- Дюран показал также, что разложение траектории утки 
по =-тени само по себе имеет тип Жевре | и этим вновь доказал существо- 
вание уток (с помощью неполного преобразования Лапласа). Его теория 
приложима к семейству более общих систем, чем уравнение Ван дер Поля. 
Это семейство содержит такие классические системы, как системы типа 
брюсселятора. 


3.4. 4-разностные уравнения 


Алгебраическим линейным уравнением в (конечных) разностях в ком- 
плексном поле называется функциональное уравнение вида 


An (a) f(a + п)... + а1(2) fla + 1) + ao(x) f(x) = 0, 


где а; — многочлены (с комплексными коэффициентами), a f — неизвестная 
функция. Например, 

2+1 -=Л (2) =0- 
разностное уравнение, допускающее в качестве решения функцию 
f(a) = Г(т). 

От дифференциального уравнения к разностному можно перейти, за- 
меняя инфинитезимальный автоморфизм 4/4х на автоморфизм переноса 
т->т-1. Это гомография римановой сферы, допускаюшая со в каче- 
стве неподвижной точки. Можно использовать и более «общую» гомогра- 
фию, допускающую две различные неподвижные точки, например, x +> gx 


64 РАСХОДЯЩИЕСЯ РЯДЫ И АСИМПТОТИЧЕСКИЕ ТЕОРИИ 


с ненулевым комплексным 4 и неподвижными точками 0 и со. Тогда при- 
ходим к понятию 4-разностного (линейного алгебраического) уравнения: 


dn (x) f(q"a) +... + ал (т) (ат) + ao(a) f(x) = 0. 


Можно доказать, что формальные ряды, являющиеся решениями ана- 
литических разностных уравнений, даже нелинейных, имеют тип Жевре: 
это аналог теоремы Мелле [GL]. Напротив, для 4-разностных уравнений 
ситуация совсем иная. В этом случае, вообще говоря, оценок Жевре не бу- 
дет, будут только оценки 4-Жевре (по терминологии, недавно введенной 
Ж.-П. Безивином (J.-P. Bezivin)). 

Мы ограничимся тем, что происходит в случае |9| 52 1. Имея возмож- 
ность заменять 4 на q~/, мы можем предполагать, что |4| > 1. 


ре +00 
Определение. Пусть f(x) = >’ anx” — формальный ряд. Говорят, 
п=0 
что f имеет 4-тип Жевре вида 3, если для подходящих С, A, s > 0 


sn(n+1) 


lan] < СЧ ? А". 


Еще говорят, что вид 3 оптимален, если не существует оценки такого же 
рода для s’ < 5. 


С каждым Qg-pa3HOCTHbIM алгебраическим линейным оператором ТГ 
можно связать многоугольник Ньютона М№(Т), аналогично случаю диффе- 
ренциальных операторов. Имеем следующий результат. 


Теорема (Ж.-П. Безивин). //усть 


Г(Л) = а. ( =) (ат) +... + a1(x) (ах) + ao(x) f(z) =0 (1) 


4-разностное линейное алгебраическое уравнение, где |4|>1. Тогда 
существует конечное число положительных действительных чисел 
31 <... < Sp, заданных сторонами многоугольника Ньютона N(T) и та- 
ких, что каждый формальный ряд a являющийся решением (1), обладает 
следующим свойством: f или сходится, или имеет q-mun Жевре оптималь- 
ного вида с одним из S;. 


Существуют 4-разностные алгебраические линейные уравнения второ- 
го порядка с решениями в виде расходящихся формальных рядов (Адамс). 
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Так что существуют расходящиеся решения 4-разностных уравнений, не 
допускающие оценок Жевре. 

Вот простой пример. Определим оператор og условием oy f(x) = 
= 7 (ах). Обозначим 


a + со n(n+1) 
CeO а 2 a 
п=0 
п(п- 1) (n+1)(n+2) 
Тогда zog(q 2 a™t!)=q 2 et?) откуда 
(хо + 1)Q(z, q) = o Q(x, q) + Q(z, Gq) =a. (2) 


Оператор то. — аналог дифференциального оператора х?4/4х; уравне- 
ние (1) — д-аналог уравнения Эйлера 


уу +у=х 
~ +00 
a ряд Q(x, q) — 9-аналог ряда Эйлера S* (—-1)"nla"t!. При gq <1 
n=0 


ряд О (т, 4) связан с 91-функцией Якоби. Из условия (og — 4)(хоа +1) = 
= 470% — то + 9 — 1 выводим, что ряд (1, 4) — решение 4-дифференци- 
ального уравнения второго порядка 


(4202 =. +а-—1)7@) =921(а*т) — af (gz) + (@—1)f (2) =0. 


С точки зрения суммируемости, случай разностных уравнений сложен: фор- 
мальные решения не будут в общем случае мультисуммируемы [E3]. Для 
а-разностных уравнений можно развивать теорию, заменяя асимптотиче- 
ские разложения Жевре на асимптотические разложения 4-Жевре, при этом 
экспоненциальное убывание заменяется на оценки вида 


f(a) < oie) | 


Имеются также 4-аналоги преобразований Бореля и Лапласа и К-суммиру- 


емости: числа n! — это моменты функции ©": 
+ со 


Г(®в 1) =и! = ] uve “du. 
0 
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v(v+1) 
Вычислим моменты функциид 2 для действительных g > 1: 


Te v(v+1) 
Ln = few . du, 
0 


log u 


и и = q’. Находим 
log q 


со = 


_1 n(n+1) 
ln = М 2710544 84а ? 


Заменяя u на &/т, получаем 


со 
Е 
Tinto" = ма" = ] &Пе dE = L(E”)(z), 
0 
где С — преобразование Лапласа и 


1 n(n+1) Cee w(w+1) 
ina“ ово 80. = fer ae 
0 


_ log€ —logz log €/ax 


сш = . Имеем 
log q log q 
nD ОВ +°° ay 41) 
2 = —____. по 2 dé. 
4 АЕ / 9 - 
Так что естественно будет определить 9-преобразование Лапласа как 
(a Te w(w+l) 
Со) = —<=€ ] а 2 d 
gh G(x) ЕЕ ) Ф(Е) - 
сш = ВЕ Отсюда 
log 4 
n(n+1) 


р) nm+1 __ т, 

q unr’ = gh (E")(2). 

Все эти задачи изучаются в настоящее время. (Я предполагаю также, что 
существует 4-вариант теоремы Мелле в нелинейном случае.) 
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3.5. Множественность естественных процессов суммирования, ветви 
функций и последнее письмо Эвариста Галуа 


Как можно было видеть по предыдущему изложению, возникает BO3- 
можность приписать расходящемуся ряду «естественную» сумму и ради- 
кально улучшить асимптотическую теорию Пуанкаре, заменив ее точной 
теорией. Так что выполняется программа, намеченная DY. Борелем [Во2 |: 


Расходящиеся ряды могут сослужить большую службу не только 
с формальной точки зрения (в чем никто никогда не сомневался) и с точки 
зрения приближенных вычислений (асимптотические ряды), но в некото- 
рых случаях их еще и можно точно вычислить. Расходящийся числовой ряд 
может иметь определенное значение. 


Но при более внимательном рассмотрении ситуация кажется менее 
идиллической. Мы видели, что в действительности в некоторых случаях по- 
является не одна, а несколько различных естественных сумм. Д. Дюмон во 
введении к своей книге [Ри] цитирует приведенный выше текст и характе- 
ризует мнение Э. Бореля как оптимистичное в сравнении с Г.Х. Харди [H1]: 


Разные методы могут давать разные суммы для одних и тех же 
рядов... 


С моей точки зрения, феномен множественности естественных сумм 
это: 


® менее удивительно, чем если бы этого не было; 


® больше заметных преимуществ, чем неудобств. 


Э.Борель? и Г.Х. Харди” очень хорошо понимали то, что, по мо- 
ему мнению, является одной из фундаментальных причин этого феномена: 
множественность аналитических продолжений. Больше всего оба они на- 
стаивали на том факте, что хорошая теория суммирования должна быть 
основана на изучении аналитических продолжений. 


23Здесь важно сделать существенное замечание: если аналитическая функция Ф(=) не од- 
нозначна, то предыдущая теория ведет к тому, что с расходящимся рядом Ф(20) связывает 
несколько разных значений, даже бесконечных [Bol, chap. 4, 6.3, р. 153]. 

2 Если У`аих" сходится при малых т и определяет функцию f(x) комплексного пере- 
менного х, однозначную и регулярную в открытой и связной области, содержащей начало 
координат и точку © = 1, и если f(x) = 3, тогда мы назовем з ©-суммой ряда >) ап. Значе- 
ние $ может естественным образом зависеть от выбора области. 
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Г.Х. Харди называет G-MeTOJOM метод суммирования с помощью ана- 
литического продолжения, на который мы ссылались в 1.5: 


...тогда мы называем ©-суммой сумму ряда У ` ат. Значение з может 
естественным образом зависеть от выбора области. 


Он также писал по поводу идей Эйлера о суммировании, о которых мы 
говорили выше, в 1.1: 


Невозможно точно установить принцип Эйлера, не имея ясного пред- 
ставления о функциях комплексного переменного и об аналитическом про- 
должении. 


Этот фундаментальный феномен появился уже при изучении аналити- 
ческого продолжения сходящегося ряда вне его области сходимости. Чтобы 
увидеть это, вернемся к уже изученному в 1.6 примеру: 


(1-0 ,.:0-00- 
р 2\2 ae 2 
f@)=T+ 524 4 


ae 
При |x| < 1 сумма ряда f равна f(z) = УТ- т. Если бы надо было 
суммировать ряд 


то можно было бы одинаково использовать аналитическое продолжение 
вдоль луча, исходящего из начала координат с аргументом 7 — = (с малым 
5 > 0) и соединенного с дугой, «убывающей» до —2, и аналитическое про- 
должение вдоль луча, исходящего из начала координат с аргументом 7 + = 
и соединенного с дугой, «возрастающей» до —2: это и есть выбор ветви 
или двойственность. В первом случае получаем 7, во втором — —$. Сравне- 
нию двух методов суммирования отвечает, в случае ряда р преобразование 
f => —/ (действие монодромии вблизи особой точки —1, порождённое осо- 
бым направлением IR” с аргументом 7). 

Феномен Стокса, открытый Стоксом при изучении уравнения Эй- 
ри (1.4), полностью аналогичен только что описанному нами феномену 
смены ветви (как мы отметили выше, эту аналогию можно найти уже в ме- 
муаре Croxca)*>. С этой точки зрения, феномен Стокса можно описать так. 


2>Расходящиеся ряды обычно делят на два класса, в зависимости от того, все ли время 
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ем + со 
Пусть f(z) = >> anx” — формальный ряд, суммируемый (или, в б0- 
n=0 


лее общем случае, мультисуммируемый) по любому направлению, близко- 
му к направлению а, за исключением (особого) направления с. Применяя 
к ряду - операторы бокового суммирования 5. и 51, мы получаем две 
разные суммы ]/ и ft. Здесь также имеет место выбор ветви, то есть 
двойственность. Работая с подходящей дифференциальной алгеброй, опре- 


делим автоморфизм Стокса, ассоциированный с особым направлением о: 


Ве 


В подходящем формализме (см. [MR2] и [MR3]) этот автоморфизм диффе- 
ренциальной алгебры интерпретируется как монодромия вокруг бесконечно 
близкой к началу координат особой точки, определённой вдоль особого на- 
правления с. 

Описание феномена Стокса, данное Стоксом для формальных реше- 
ний уравнения Эйри на бесконечности, очень похоже на только что рас- 
смотренное нами. Грубо говоря, он действует следующим образом: строит 
базу В решений уравнения Эйри, используя сходящиеся разложения в нача- 
ле координат (восходящие ряды). Затем он вновь суммирует формальные 
решения на бесконечности NO разным направлениям (нисходящие ряды) 
с помощью экспоненциально точного численного метода (вариант сумми- 
рования до наименьшего члена). Этот метод достаточно точен, чтобы позво- 
лить Стоксу выразить суммы формальных рядов на бесконечности в базе В 
с помощью произвольных констант; сегодня это выражение известно как 
формула связи. Тогда возникает вот какая проблема: разложения в нача- 
ле координат — целые функции (начало координат — регулярная точка), 
и поэтому монодромия вокруг начала координат тривиальна, тогда как рас- 
ходящиеся разложения на бесконечности содержат \/т и их монодромия 
(формальная монодромия) нетривиальна, что явно противоречит только что 
проведенному анализу. Именно это противоречие вызывало у Стокса рас- 
терянность в течение долгих лет, прежде чем он нашел ключ к загадке 
(см. процитированное выше письмо к невесте): 


положительны их члены или попеременно положительны и отрицательны. .., ряды предыду- 
щего вида появляются как особый вариант общего случая обращения с расходящимися рядами 
согласно степеням мнимой переменной, как неопределенные формы, при прохождении через 
которые аналитическое выражение терпит разрыв, аналогичный смене знака в радикале [Sto2, 
р. 78]. 
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... NOCKONbKY нисходящие ряды содержат радикалы, которые не появ- 
ляются в восходящих рядах, мы можем видеть а priori, что произвольные 
константы должны быть разрывны. 


Описанная нами двойственность появляется также как разрывность 
констант связи ири пересечении особой прямой: для особой прямой точ- 
ность численного метода суммирования Стокса недостаточна для вычисле- 
ния произвольных констант (одна из них теряется; в действительности речь 
идет о точном вычислении экспоненциально разбегающегося решения, а для 
этого надо располагать достаточной экспоненциальной точностью; вдоль 
особой прямой разбегающееся решение слишком мало, чтобы рассматри- 
вать его численно с помощью используемого метода!). Продолжим анализ 
Дингла [Di, chap. 1, p. 7]: 

Лучи Cmoxca для асимптотических рядов определяются теми фазами, 
для которых ряд (включая его множитель) достигает максимума экспо- 
ненциального преобладания над связанной с ним функцией. 


(Разрывность произвольных констант появляется в момент максималь- 
ного преобладания экспоненциально доминантного символа над экспонен- 
циально рецессивным.) 

Только что развитая нами точка зрения заметно отличается от тради- 
ционного подхода к феномену Стокса”°. Согласно этому подходу феномен 
состоит в смене преобладания среди двух экспонент; итак, его можно на- 
блюдать вдоль линий осцилляции или линий Стокса. Напротив, в нашем 
описании (и у самого Стокса) феномен проявляется на особых прямых 
(которые иногда называют линиями анти-Стокса). Источник этого разли- 
чия — в противопоставлении традиционного видения расходящихся рядов 
как асимптотических и концепции расходящихся рядов как кодирующих 
точные решения. В одном случае акцент делается на асимптотике в смыс- 
ле Пуанкаре (и фундаментальная природа феномена не воспринимается), 
в другом случае используется точная асимптотика (см. также [СМР]. 

Продолжим сравнивать смену алгебраических ветвей и смену ветвей 
в феномене Стокса. Если между этими двумя феноменами есть глубокая 
аналогия (речь идет о двух случаях преобразовании Галуа: об автоморфиз- 
мах дифференциальных алгебр), то есть также и радикальные различия. 


Только Дингл, по-видимому, устраняется от этого подхода и продолжает оригинальные 
идеи Стокса. Впрочем, он называет линиями Стокса то, что другие авторы называют линиями 
анти-Стокса (у нас: особые прямые). Это соответствует смыслу одной из его центральных 
идей: поиску точной асимптотической теории. Его полные асимитотические разложения — 
прообраз трансасимптотических разложений. 
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Например, матрица смены ветви для (1+ 2, V1 + 23) равна 


©“). 


тогда как такая же матрица феномена Стокса St, для уравнения Эйлера 
(соответственно для уравнения Эйри) имеет вид 


(0 1): 


где В — ненулевое комплексное число. 

Вторая матрица унипотентна, так как феномен Стокса не обнаружи- 
вается асимптотически при переходе через особую прямую: после перехода 
асимптотические разложения решений не изменились (что совершенно не 
так в случае смены алгебраических ветвей). Можно показать — в основном, 
с помощью тех же рассуждений, — что операторы Стокса всегда унипо- 
тентны. Отсюда фундаментальное следствие: операторы Стокса St „ — это 
инфинитезимальные порождающие (логарифмы), обозначаемые До: 


= eda, 


где оператор Nir ot производная Галуа, то есть производная, коммутиру- 
ющая с обычной, именуемая ориентированной посторонней производной. 
Это позволяет изучать феномен Стокса инфинитезимально. Эту точку зре- 
ния открыл и систематически изучал Жан Экалль: Восходящие функции, По- 
стороннее дифференциальное исчисление, Ускорение [Е1], [E2], [E3], [E4], 
[Са] и [CNP]. 

Накануне дуэли, на которой ему предстояло умереть, Эварист Галуа 
писал в своем последнем письме, адресованном своему другу Огюсту Ше- 
валье [G1]: 


Знаешь, мой дорогой Огюст, эти темы — не единственные, которые 
я исследовал. С некоторых пор главные мои размышления были направлены 
на приложение теории двойственности к трансцендентному анализу. Хо- 
телось бы a priori увидеть связь между трансцендентными величинами, 
или трансцендентными функциями, увидеть, какие замены можно было 
бы сделать, какие величины можно было бы подставлять вместо данных 
так, чтобы эта связь не могла исчезнуть. Это заставляет сразу признать 
невозможными многие выражения, которые можно было бы искать. Но 


72 РАСХОДЯЩИЕСЯ РЯДЫ И АСИМПТОТИЧЕСКИЕ ТЕОРИИ 


у меня нет времени, а мои идеи на этой, вообще говоря, обширной почве 
еще не очень развиты. 


Значение этого текста, по-видимому, долгое время оставалось довольно 
загадочным. Но сегодня мне кажется возможным прояснить большую его 
часть. В своем предисловии к /1олному собранию сочинении Галуа [G3] Жан 
Дьедонне писал: 


...HO остается место для раздумий о том, что он должен был быть 
очень близок к идее «римановой поверхности», связанной с алгебраической 
функцией, и что такая идея должна была стать фундаментальной в его 
исследованиях, которые он называл «теорией двойственности». ... 


Эта гипотеза несколько проясняет вопрос, но если внимательно перечи- 
тать слова Галуа, то кажется очевидным, что идея римановых поверхностей 
в лучшем случае касается только одной части его «теории двойственности». 
Я думаю, что Софус Ли, Эмиль Пикар и особенно Жюль Драш угадали, в ка- 
ком направлении в действительности шли последние исследования Галуа. 
В своем предисловии к Полному собранию сочинений Галуа [G2] Эмиль 
Пикар писал: 


..OH ...создал бы, в ее существенной части, теорию алгебраических 
функций одной переменной в том виде, в каком мы ее знаем сегодня. Раз- 
мышления Галуа шли еще дальше; он заканчивает свое письмо словами 
о приложении теории двойственности к трансцендентному анализу. Мож- 
но почти угадать то, что он под этим понимал, и на этой обширной, как 
он говорил, почве еще и сегодня остается много открытий, которые пред- 
стоит совершить... 


Очевидно, нам остается угадать, что хотел сказать Пикар! Если вспо- 
мнить, что он один из основателей дифференциальной теории Галуа (тео- 
рия Пикара-Вессио), то кажется достаточно разумным думать, что, по его 
убеждению, у Галуа были некоторые идеи относительно этой теории. Во 
всяком случае, в этом был убежден Жюль Драш, написавший в конце своей 
диссертации: 


Мы были бы счастливы, если бы наша работа могла привлечь внимание 
к нескольким строкам, которыми заканчивается письмо Галуа, и если бы ее 
можно было рассматривать как первую попытку объяснить выраженные 
там мысли. 

Вместе с тем, несколькими годами раньше, в 1895, во время своей 
лекции об этой работе Галуа в честь столетия Высшей Нормальной Школы 
Софус Ли выступил точно в том же духе [L]: 
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Последние указания в научном завещании Галуа не менее замечатель- 
ны, хотя он выразил их в такой смутной форме, что с трудом можно 
угадать его мысль. Деиствительно, нельзя сомневаться, что Галуа наме- 
ревался исследовать не только группы подстановок, но и группы преобразо- 
ваний с совершенно общей точки зрения и что он мечтал об их применении 
в анализе. Сейчас практически невозможно понять, думал он о непрерыв- 
ных группах, или о разрывных, или о тех и других одновременно; и мы 
не можем сказать, какую пользу он собирался из них извлечь, но, конечно 
же, Галуа предчувствовал значение, которое могли приобрести в других 
областях те идеи, которые привели его к столь блистательному успеху 
в теории алгебраических уравнений. 

Со своей стороны, вот уже несколько лет как я читаю письмо Галуа, 
и мне сразу показалось очевидным”’, что идеи Галуа по теории двойствен- 
ности предварили дифференциальную теорию Галуа и даже — в некоторой 
мере — то, что я недавно читал по этой теории. Невозможно узнать, имел 
ли Галуа понятие о феномене Стокса и о его природе, связанной с теорией 
Галуа, что я счел очевидным в [Ка3]. (Изучение фрагментов вычислений, 
найденных в его бумагах [G3], не позволяет делать такие выводы.) Напро- 
тив, я, каки Жюль Драш, уверен, что Эварист Галуа понимал, что некоторые 
преобразования в комплексном анализе (двойственности. .. ) имеют «приро- 


ду Галуа», как, например, «рекалибровка» экспонент (ее прототипом будет 
1 1 


замена е? на Ле? с фиксированным ненулевым комплексным числом Л BO 
всех формулах). Очевидно, эта уверенность очень субъективна, и я не в со- 
стоянии обосновать это утверждение, что не мешает мне верить в него... Во 
всяком случае, я считаю, что развитие исследований в этой области должно 
в конце концов сделать «очевидным» понимание последних размышлений 
Галуа. 

Для дифференциальных уравнений дифференциальная теория Галуа 
играет ту же роль, что и классическая теория Галуа для алгебраических 
уравнений. Чтобы как следует понять этот предмет, лучше всего обратиться 
к оригинальной идее Галуа, выраженной им в особенно ясной манере [G2], 
[G3]: 


27Во время этого чтения я не знал ни о тексте Драша, случайно открытом мной несколькими 
неделями позже, когда я искал ссылку по просьбе H. Камрана, ни о тексте Пикара, о котором 
мне сообщил Д.Беннекин после обзора моей интерпретации последнего письма Галуа на 
моем семинаре в Страсбурге, ни о тексте С.Ли, о котором мне сообщил Ж. Б. Бост после 
опубликования первой версии этого текста. .. 
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Предложение Г 3b. Пусть дано уравнение, имеющее т корней 
a, 6, с,... Всегда найдется группа перестановок букв а, 6, Cc, ..., удовле- 
творяющая следующему условию: 

1) пусть каждая функция корней, инвариантная при подстановках из 
этой группы, будет рационально-известной; 

2) обратно, пусть каждая рационально определенная функция корней 
будет инвариантна при этих подстановках. 


Дифференциальная теория Галуа была открыта Пикаром и Вессио [Pi], 
[У]. Заинтересованный читатель может просмотреть замечательное вве- 
дение [Ка]. Так же можно использовать [MR2], [Каб] и более техни- 
ческие статьи [Ra4], [Mi], чтобы узнать больше о связи расходящих- 
ся рядов с дифференциальной теорией Галуа. Здесь мы удовлетворим- 
ся несколькими указаниями и при этом ограничимся линейным случа- 
ем. Пусть (К, д) — дифференциальное тело нулевой характеристики, 
а С = {c € K; дс = 0} — тело его констант. (Например, для (C(x), а/ах), 
(C{x}, х?а/ат) и (C{{x}}, х?а/ат) телом констант будет С.) 

Пусть D = a,(d/dx)" +...+a,d/dz+ag — дифференциальный опера- 
тор с коэффициентами из А. Связанное с D расширение Пикара-Вессио L — 
дифференциальное надтело для А, дифференциально порожденное над К 
фундаментальной системой решений D; при этом тело констант С' остается 
неизменным. Если С’ алгебраически замкнуто, то такое расширение су- 
ществует и единственно (с точностью до изоморфизма). По определению, 
дифференциальная группа Галуа оператора D над телом К — это группа 
Gal к(р) = Autx(L) К-автоморфизмов дифференциального тела Г. Эле- 
мент о из Gal к (D) — это векторное С-пространство решений D в Г. Отсю- 
да выводим (фиксируя базис этого пространства) представление Gal x (D) 
в GL(n; С). Полученная подгруппа группы GL(n; С’) — алгебраическая (то 
есть она определена многочленами п? переменных). 

Дифференциальные группы Галуа можно вычислять, зная монодромию 
(формальную или нет), упомянутую выше «рекалибровку» экспонент и фе- 
номен Стокса. 

Чтобы закончить это направление, процитируем довольно удивитель- 
ный результат, полученный К. Мичи (С. Mitschi) с помощью суммирования 
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расходящихся рядов [МИ?"°. Рассмотрим функцию 


+ со 1 
_= 7 
А (= ] ee ae 
0 
Sie, ae 

Эта функция — преобразование Лапласа для т 2e~* , она удовлетворя- 
ет линейному алгебраическому дифференциальному уравнению седьмого 
порядка: 
1 
2 


Если в этом дифференциальном уравнении произвести замену переменной 


Dek А OO ik Bek =o. 


Е (41) (ветвление), то полученная из А функция U удовлетворяет 


обобщенному конфлюэнтному гипергеометрическому дифференциальному 
уравнению 


Вл 0 = 0, 
где D714 = 2(O+ a Е (a ~ tr) ид = zd/dz. 


Можно показать, что для D7 1 дифференциальная группа Галуа равна?” 


Gal C(x)(D7,1) = Go x 0/24. 


Монодромия оператора Dx очевидно тривиальна; отсюда: дифферен- 
циальная группа Галуа для к равна 


Gal o(2) DK = Go. 


Связь между «особой функцией» А и группой исключений Со по- 
хожа на связь, существующую между функциями Эйри А или В и груп- 
пой SL(2; С) [MR2]. 

В конце этой работы мне остается отметить, что я не смог рассказать 
обо всех приложениях расходящихся рядов, даже ограничиваясь случаем, 


28 Несколько раньше H. Кац [Kat2], [Kat3] получил тот же результат с алгебраической точки 
зрения. 

29 Группа G'2 (или скорее ее представление в размерности 7) — это алгебраическая подгруппа 
комплексной размерности 14 специальной ортогональной группы SO(7; С). Действительно- 
значная разновидность этой группы обладает хорошими свойствами: это группа автоморфиз- 
мов алгебры октав Кэлли и — топологически — слой базы сферы 56, чьи слои — это слои 
базы 55 и слоя 93 [Pos]. 
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когда имеют место оценки Жевре и/или экспоненциально малые поправ- 
ки. (Среди наиболее значительных пропусков я бы процитировал недавние 
замечательные работы о стационарном уравнении Шредингера в полуклас- 
сической «точной» постановке [Vo], [CNP], как и принадлежащее Шеффке 
и Фолькмеру доказательство несуществования изокордных овалов с двумя 
«центрами».) Другие области математики или физики 


® адиабатические инварианты, 


® квантовая теория поля (инстантоны, ренормалоны, ряды возмущений 
И Т. Д.), 


® «Ттермодинамическое» изучение задачи о коммивояжере, 


® СОЛИТОНЫ И Т. Д., 


очевидно, зависят от этой последней проблематики, хотя на данный мо- 
мент здесь немного точных теоретических результатов. Остается огромное 
количество работы! 
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